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SECTION  I 
INTRODUCTION 

The  desirable  properties  of  magnetic  bubble  memories  for  Air  Force 
remote-vehicle  application  are  well-known.  The  properties  of  small 
physical  size,  low  power,  light  weight,  radiation  tolerance,  and  low  cost 
are  combined  with  nonvolatility  and  block-access  capability  which  make 
magnetic  bubbles  an  ideal  substitute  for  discs,  drums,  and  modest-capacity 
tape  recorders.  Magnetic  bubble  memories  fit  precisely  into  the  computer 
hierarchy  in  place  of  presently-used  electromechanical  technologies,  but 
with  improved  operating  characteristics. 

The  major  challenge  to  magnetic  bubble  memories  has  been  to  provide 

enough  bit  capacity  in  a  small  volume  for  Air  Force  remote  vehicle 

applications.  The  first  Air  Force  Avionics  Laboratory  program  to  construct 

bubble  memories  uses  four  to  six  micrometer  diameter  bubbles  and  approxi- 
5  1 

mately  10  -bit  chips.  To  keep  the  reliability  and  lifetime  high,  only  a 
few  thousand  chips  can  be  wired  together.  Thus  the  practical  limit  to 

g 

memory  capacity  using  four  to  six  micrometer  bubbles  is  about  3x10  bits. 

The  volume  required  for  a  bubble  memory  with  this  high  a  capacity  becomes 

rather  large  (several  cubic  feet),  so  the  typical  disc/drum  capacity  of 

about  1.5  x  10^  bits,  which  can  be  accomplished  in  about  one  cubic  foot, 

is  a  practical  goal  for  the  first  (large  bubble)  memory  systems.  It  is 

typical  of  emerging  technologies  that  users  are  not  satisfied  with  easily 

attainable  goals,  and  magnetic  bubbles  are  no  exception.  Even  before  the 

work  to  construct  large  bubble  memories  was  started,  the  need  for  lower 

volume  and  higher  total  bit  capacity  was  apparent.  In  order  to  attain  a 

g 

capacity  of  2  x  10  bits  in  no  more  than  0.5  cubic  feet,  a  major  advance 
in  magnetic  bubble  technology  was  needed.  Reliability  considerations 
continue  to  dictate  that  no  more  than  a  few  thousand  chips  be  used,  and 
system  volume  considerations  dictate  that  the  number  of  coil  sets  be  kept 
as  low  as  possible.  Thus  the  chip  capacity  must  be  the  order  of  10^  bits. 
While  four  to  six  micrometer  diameter  bubbles  could  be  used  in  a  chip  of 
this  capacity,  the  chip  would  be  several  square  inches. 
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Maintaining  a  uniform  bias  and  rotating  magnetic  field  over  this 

area  is  a  complex  problem.  A  simpler  approach  is  to  use  small  (two 

micrometer  diameter)  magnetic  bubbles.  This  is  the  second  group  of  AFAL 

program  aimed  at  million-bit  or  larger  chips  and  systems  with  capabilities 

2 

up  to  a  few  billions  of  bits.  Such  memory  systems  could  fill  the  position 

in  computer  hierarchies  now  held  by  moving-head  discs  and  also  many 

recorder  application.  Large  improvements  in  performance,  ruggedness, 

reliability,  and  lifetime  should  be  achieved  by  utilizing  small  bubble 

technology.  Recording  and  processing  sensor  data,  as  from  side-looking 

radar,  forward-looking  infrared,  or  visual  imagery,  can  require  much 

larger  data  memories.  Section  II  presents  the  analysis  of  the  in/out 

bit  rate  and  total  on-board  storage  capacity  that  would  be  needed  for  a 

specific  (but  hypothetical)  border  surveillance  task.  It  is  shown  that 
1 2 

10  bits  storage,  with  an  in/out  data  rate  of  50  Mbit/second  is  needed 
on-board  a  remote  airplane.  This  capacity  would  require  over  a  hundred 
cubic  feet  volume  if  implemented  using  small  bubbles,  and  it  is  evident 
that  another  major  advance  in  magnetic  bubble  technology  is  needed.  The 
approach  chosen  is  to  use  lattice  file,  whole  water  chips  which  could 

O 

achieve  at  least  2x10  bits  per  chip  and  780  Kbit/second  in/out  rates. 
Multiplexing  sixty- four  of  these  chips  would  give  the  desired  50Mbit/ 
second  data  rate.  This  research  is  the  subject  of  Contract  No.  F3361 5-76- 
C-1198. 

The  in-house  portion  of  this  basic  research  has  three  aims:  (1)  to 
perform  the  system  concepts  analysis  of  what  the  Air  Force  will  need  to 
perform  a  mission  in  the  mid-1980's  time  frame,  (2)  to  ensure  that 
contractor  efforts  actually  contribute  to  achievement  of  the  goals 
determined  during  the  system  concepts  analysis,  (3)  to  contribute  new 
ideas  to  the  contractor  and,  where  needed,  to  help  the  contractor  overcome 
impasses  in  his  approach.  Section  III  addresses  a  specific  impasse: 
optical  methods  are  inadequate  to  characterize  magnetic  bubble  materials 
when  the  thickness  and  bubble  diameter  are  only  three  or  four  times  the 
wavelength  of  visible  light.  A  new  non-optical  technique  for  characterizing 
magnetic  bubble  materials  is  described  in  this  section  and  compared  to 
the  more  traditional  bubble  materials  characterization  techniques.  Finally, 
the  computer  program  used  in  performing  this  new  analysis  is  listed  in  the 
Appendix  with  instructions  for  use. 
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SECTION  II 

GEOMETRICAL  ANALYSIS  OF  AN  AIR  FORCE  TACTICAL  MISSION  SCENARIO 

1.  INTRODUCTION 

Military  power  sometimes  has  been  equated  with  numbers  of  planes, 
tanks,  ships,  guns,  and  men,  but  there  is  one  critical  element  in  addition  - 
that  is  information.  All  weapons  are  completely  useless  without  specific 
information  on  where  the  targets  are  located.  It  is  the  function  of 
reconnaissance  to  explore  the  enemy  territory  and  gain  information  about 
the  location  of  possible  threats  or  targets.  This  information  may  be 
obtained  in  the  form  of  photographs,  radar  imagery,  infrared  imagery,  or 
other  sensory  data.  Modern  electronics  usually  reduces  this  information 
into  digital  form  to  facilitate  handling  by  digital  computer.  The 
transmission  of  this  information  by  various  types  of  communication  systems 
and  the  manipulation  of  the  information  into  a  form  usable  for  command 
and  control  decision  making  is  the  function  of  data  processing.  Regardless 
of  how  the  information  originated,  it  ultimately  will  reside  at  some  state 
of  processing  in  a  computer  mass  memory.  It  is  the  purpose  of  this  analysis 
to  show  how  the  development  of  smaller,  lighter,  higher  capacity,  faster 
mass  memories  has  a  direct  and  immediate  effect  on  the  use  of  tactical 
weapons.  Just  as  information  is  a  critical  ingredient  of  military  power, 
so  are  the  memories  which  contain  this  information. 

2.  GEOMETRICAL  SCENARIO  ANALYSIS 

Military  applications  for  solid-state  mass  memories  include  such  areas 

as:  radar  video  recording,  electronic  intelligence,  reconnaissance, 

surveillance,  and  general  tape  recorder  replacement.  For  purposes  of 

illustration  consider  one  of  the  critical  tactical  problems  discussed 

by  Dr.  Mai  comb  Currie,  director  of  DDR&E ,  in  his  statement  before  congress 

3 

on  the  Fiscal  1977  budget  (Figure  1).  That  problem  is  how  to  counter  a 
"blitzkrieg"  of  massed,  highly  mobile,  heavy  armor  such  as  we  may  some  day 
face  ’n  Europe.  The  solution  to  this  problem  was  suggested  by  Dr.  Currie 
through  the  use  of  improved  command/control  "force  level  multipliers" 
together  with  integrated  battlefield  surveillance,  target  acquisition. 
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high  altitude  sensors,  stand-off  synthetic  aperture  radars,  secure  com¬ 
munications,  computers  capable  of  real  time  processing  of  sensor  data, 
and  large  numbers  of  surgically  accurate  weapons  with  less  sophistication 
and  low  cost. 

Small  solid-state  mass  memories  are  one  of  the  key  components  along 

with  high-speed  computers  and  low-cost  strike  Remotely  Piloted  Vehicle 

(RPV's)  to  the  implementation  of  this  strategy.  In  order  to  store  data 

on  an  area  representing  the  size  of  West  Germany  (2.6  x  10  sq.  ft.), 

for  example,  with  a  resolution  of  1  yard  and  1  bit  of  gray  scale  would 
1 2 

require  0.29  x  10  bits  of  memory.  All  present-day,  on-line  memories 
with  this  capacity  require  critical  mechanical  alignments  (and  are  thus 
fragile).  They  typically  occupy  at  least  a  thousand  cubic  feet  of  air- 
conditioned,  vibration-free  computer  room  space.  Within  these  constraints 
it  is  a  remarkable  achievement  that  three  versions  have  found  use.  If  one 
assumes  a  hypothetical  number  of  10,000  targets  in  this  area,  free  to  move 
at  up  to  50  mph,  and  that  we  wish  to  know  their  position  accurately  to 
within  100  feet,  we  must  update  the  data  on  all  targets  every  1.36  seconds 
and  the  memory  must  be  capable  of  an  access  time  no  greater  than  1.36  x  10 
seconds.  To  fit  aboard  a  modern  fighter  aircraft  would  require  that  the 
memory  be  no  more  than  a  couple  of  cubic  feet  in  size  and  weigh  less  than 
a  few  hundred  pounds. 

3.  CURRENTLY  AVAILABLE  TECHNOLOGIES 

A  typical  capacity  for  present-day  large  data  base  mass  memories  is 

12  11 
about  10  bits.  A  large  reel  of  magnetic  video  tape  will  hold  about  10 

bits  and  some  means  of  mechanically  changing  reels  must  be  provided  to 
1 2 

reach  10  bits.  All  present-day,  on-line  memories  with  this  capacity 

require  critical  mechanical  alignments  (and  are  thus  fragile).  They 

typically  occupy  at  least  a  thousand  cubic  feet  of  air-conditioned, 

vibration- free  computer  room  space.  Within  these  constraints,  it  is  a 

remarkable  achievement  that  three  versions  have  found  use.  The  charac- 

1 2 

teristics  of  currently  available  10  bit  mass  memories  are  summarized 
in  Tabl e  1 . 
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TABLE  1 

CURRENTLY  AVAILABLE  1012  BIT  MASS  MEMORIES 


Precision  Instruments  Ampex  Terabit 
12  12 

System  Capacity  1  x  10  bits  1.5  x  10  bits 

6  5 

Data  Density  3x10  bits/sq.  cm.  1  x  10  bits/sq.  cm. 

5  7 

Recording  Area  4  x  '0  sq.  cm.  4  x  10  sq.  cm. 
Access  Time  5  seconds  20  seconds 


IBM  3850 
2  x  1012  bits 

4 

5  x  10  bits/sq.  cm. 

4  x  10^  sq.  cm. 

5  to  8  seconds 


All  of  these  systems  use  some  kind  of  plastic  tape  with  information  stored 
magnetically  or  optically  (as  laser  burned  holes)  with  mechanical  access. 
The  slow  access  times,  together  with  the  very  large  physical  size,  preclude 
their  use  in  any  practical  airborne  or  spaceborne  system  foreseen  today. 


4.  EMERGING  TECHNOLOGIES 

Prospects  for  improving  mass  memories  in  the  future  look  extremely 

1 2 

promising.  The  four  leading  candidates  for  10  bit  solid  state  memories 
in  the  next  ten  years  are:  semiconductor  Charged-Coupled  Devices  (CCD's), 
magnetic  bubbles,  optical,  and  electron  beam  accessed  memories.  Semi¬ 
conductor  CCD's  were  selected  to  illustrate  the  state  of  the  art  in 
semiconductor  memories.  The  Mnemonics  65K  CCD  chip  is  the  largest 
capacity  semiconductor  memory  chip  currently  on  the  market.  Data  on 
this  chip  together  with  projections  for  CCD  technology  presented  by 
Dr.  W.  Kosonocky  at  the  Stanford  Research  Institute  Symposium  on  Advanced 

4 

Memory  Concepts  are  listed  in  Table  2.  The  main  disadvantage  of  today's 
CCD  technology  is  the  very  large  standby  power  necessary  to  constantly 
recirculate  and  refresh  data,  since  this  is  a  volatile  type  of  memory. 
Substantial  progress  can  be  anticipated  over  the  next  ten  years  if  one 
assumes  that  one  micron  design  rules  are  common  practice,  that  five  inch 
whole  wafer  technology  is  available  with  yields  of  at  least  20%,  and 
that  some  combination  of  CCD  techniques  with  non-volatile  Netal  Nitride 
Oxide  Semiconductor  (MN0S)  technology  will  become  available  to  produce 
a  nonvolatile  high-capacity  CCD  which  would  not  require  standby  power. 

If  these  three  assumptions  can  be  realized  without  increasing  the  cell 
size  too  much,  then  CCD's  may  become  much  more  competitive  in  the  future. 
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TABLE  2 

SEMICONDUCTOR  MEMORIES 


Type 

Volatile/Non-Volatile 

Category 

Cell  Size 

Chip  Size 

Access  Time 

1/0  Data  Rate 

Volume 

Weight 

Power 

Temperature  Range 
Cost 


Current  Capabilities  Projected  Copobilities  in  1985 

CCD  (NMOS  -  Double  Polysilicon)  CCD 


Volotile 

Read/Write 

22.7  x  22.7  Microns 

218  x  235  mils 

0.4  to  2  mSec 

5  Mbps 

2404  cu  ft  for  10"  Bits 
269,248  lbs  for  10"  Bits 
246,154  Watts  Standby  for  10" 
0*  to  70*  C 
10  m$/Bit 


Volatile 
Read/Write 
4x4  Microns 
5  in  Diam  Wafer 
0.1  to  20  mSec 
10  Mbps 

14  cu.  ft  for  10" 

1570  lbs  for  10" 

1  -  10  KW  Standby  for  10" 
C*  -  70*  C 
.1  to  10  m  4/Bit 


Present  and  projected  characteristics  for  optical  memories  are  based 

primarily  on  the  Precision  Instruments  Unicon  190  system  as  described  by 

Kaczorowski6  and  Dell6  (Table  3).  This  system  utilizes  a  laser  to  burn 

1 2 

small  holes  in  a  metal  film  on  a  polyester  strip.  A  10  bit  capacity 
system  is  presently  available  at  relatively  low  cost.  However,  this 
system  requires  relatively  high  power  for  the  read  laser  and  exhibits 
very  slow  access  speeds  due  to  mechanical  motion  necessary  for  changing 
the  storage  strip.  Experimental  systems  have  been  built  without  mechanical 
motion,  but  the  costs  of  the  optics  become  prohibitive:  a  separate  read 
station  must  be  provided  for  each  storage  strip.  It  is  anticipated  that 
the  density  of  storage  in  these  systems  will  increase  by  about  a  factor 
of  ten  on  the  medium  as  one  micrometer  cell  sizes  are  achieved,  and  that 
this  will  lead  to  somewhat  faster  access  and  lower  system  volume.  Optical 
memory  systems,  however,  will  always  be  limited  by  the  laser  in  power 
efficiency  as  well  as  in  size  and  weight.  In  addition  they  require 
precision,  critical  alignments  which  are  difficult  to  maintain  in  aero¬ 
space  environments. 
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TABLE  3 

OPTICAL  MEMORIES 


Typo 

Volatile/Nonvolatile 

Category 

Cell  Size 

Chip  Size 

Chip  Capacity 

Access  Time 

I/O  Data  Bate 

Volume 

Weight 

Porrex 

Cost 


Current  Copobilltleg, 

Bit  by  Bit 
Non-Volatile 
Archival 
3.6  x  3.6  uM 
4.75  x  31.25  In 
2.8  x  10*  Bits 
10  Sec 
4  Mbps 

850  cu  ft  for  10°  Bits 
350  lbs  for  10*  Bits 
2500  V/atts  for  10"  Bits 
0.1  m<t/Bit 


Projected  Capabilities  In  1985 

Bit  by  Bit 

Non-Volatile 

Archival 

1  x  1  uM 

4  x  4  in 

10“  Bits 

1  Sec 

10  Mbps 

27  cu  ft  for  10*  Bits 
300  lbs  for  10*  Bits 
2500  Watts  for  10"  Bits 
0.06  m$/Bit 


Electron  beam  accessed  memory  capabilities  are  displayed  in  Table  4, 
where  present  capabilities  are  those  listed  by  General  Electric  Co.  and 
Microbit  Corp.  and-  1985  projected  capabilities  are  based  on  the  work  of 
G.E.  under  contract  F336T 5-76-C-T 322^ .  This  technology  has  the  potential 
of  becoming  one  of  the  least  expensive  and  fastest  of  any  of  the  current 
candidates  due  to  the  extremely  high  potential  density  of  the  data  on  the 
target  and  the  speed  with  which  an  electron  beam  can  be  deflected.  It 
does  have  the  usual  disadvantages  of  tube  technology  in  requiring  a  vacuum 
envelope  and  high  voltage  power  supplies  and  will  most  likely  require 
precision,  critical  alignments  in  the  deflection  elements.  This  tends  to 
make  it  somewhat  heavier,  more  bulky,  and  more  fragile  than  is  desired 
for  an  aerospaceborne  memory,  so  that  this  must  be  carefully  traded 
against  speed  and  cost. 


Table  5  lists  the  current  status  and  system  characteristics  of  magnetic 

bubble  mass  memories  based  on  the  Texas  Instruments  Inc.  work  under 

8 

contract  F3361 5-75-C-l 228  and  the  projected  characteristics  of  bubble- 

lattice-file  memories  based  on  work  at  Rockwell  International  under 

g 

contract  F3361 5-76-C-l 1 98  .  The  outstanding  advantages  of  magnetic 
bubble  technology  are  extremely  small  size,  weight,  and  power  with 
relatively  low  cost.  The  present  disadvantages  are  a  somewhat  limited 
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Type 
Cell  Size 
Chip  Size 
Access  Time 
I/O  Data  Rate 
Volume 
Weight 
Power 

Temperature  Range 
Total  System  Cost 


Type 
Cell  Size 
Access  Time 
I/O  Data  Rate 
System  Volume 
System  Weight 
System  Povrer 
Temperature  Range 
Cost 


TABLE  4 

ELECTRON  BEAM  MEMORY 


Current  Capabilities 
MOS  Target  (Read/ Write) 
4x4  Microns 
1  in  Wafer  for  32  Mbits 
10  —  30  uSec 
10  Mbps 

1.67  cu  ft  for  32  Mbits 
70.4  lbs  for  32  Mbits 
256  Watts  for  32  Mbits 
-55*  C  to  +125*  C 
50  m  $/Bit 


Projected  Capabilities  in  1985 
Ion  Implanted  Si  Target  (Archival)' 
0.1  uM  x  0.1  uM 
3  in  Diam  Wafer  for  10^  Bits 
30  uSec 

10  Mbps  (Single  Channel) 

3.6  cu  ft  for  1012  Bits 
700  lbs  for  1012  Bits 
256  Watts  for  1012  Bits 
-55*  C  to  +125*  C 
0.02  m$/Bit 


TABLE  5 

MAGNETIC  BUBBLE  MEMORIES 


Current  Copobilities 
Permalloy  Bor  File 
5  uM  Bubbles,.  22  uM  x  22  uM 

1.5  mSec 
2  Mbps 

2.6  cu  ft  for  100  MBit.- 
156  lbs 

51  Watts 

-25*  C  to  +125*  C 
50  m^/Bit 


Projected  Capabilities  in  1985 
Bubble  Lattice  File 
2  uM  Bubbles,  3  uM  x  3  uM 
80  uSec 
50  Mbps 

1.4  cu  ft  for  10°  Bits 
157  lbs  for  10“  Bits 
102  Watts  for  10u  Bits 
-55*  C  to  +125*  C 
0.25  m  ♦/Bit 
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temperature  range  (which  is  still  much  wider  than  mechanically-accessed 
memories  with  similar  capacities)  and  a  somewhat  slower  internal  shift 
rate  than  desired  for  some  applications  (of  course,  a  slow  internal  data 
rate  can  be  multiplexed  to  any  desired  in/out  data  rate  if  the  access 
delays  are  not  too  large).  The  projections  for  1985  are  that  many  of 
these  difficulties  will  have  been  overcome,  particularly  with  the  appli¬ 
cation  of  the  bubble  lattice  file  concept. 

1 2 

In  Table  6  the  1985  projected  characteristics  of  10  bit,  new,  mass 
memory  technologies  are  compared.  The  user's  choice  of  a  mass  memory 
technology  will  not  be  a  simple  decision.  If  one  is  primarily  concerned 
with  read/write  updateable  memories,  there  appear  to  be  only  two  choices: 
semiconductor  CCD's  and  magnetic  bubbles.  Of  these  two,  semiconductor 
CCD's  should  provide  somewhat  faster  access  times  and  should  be  capable 
of  high  in/out  data  rates  with  less  multiplexing,  but  it  appears  that 
system  volume,  weight,  and  power  will  be  much  greater  (especially  if  some 
technique  for  producing  very  large  capacity  chips  providing  nonvolatile 
data  storage  is  not  yet  available).  Magnetic  bubbles  should  be  able  to 
provide  minimum  size,  weight,  and  power  at  moderate  cost  and  may  well  be 
the  optimum  solution  to  many  user  needs.  If  archival  data  storage  (write 
once  and  read  many  times;  change  media  to  rewrite)  is  needed,  the  choice 
may  include  electron  beam  and  optical  memories.  If  high  speed  access  and 
high  data  rate  are  of  prime  importance,  the  electron  beam  memory  has  much 
to  offer.  If  low  media  cost  is  of  prime  concern  and  slow  access  time  is 
acceptabl e, then  optical  memories  may  be  best. 

5.  CONCLUSION 

We  have  accomplished  the  following  in  this  section: 

(1)  We  have  derived  a  set  of  mass  memory  specifications  from  a 
realistic  (but  hypothetical)  Air  Force  mission  scenario. 

(2)  We  have  examined  currently  available  mass  memories  to  see  if 
any  of  them  can  meet  this  specification 

(3)  We  have  projected  four  technologies  that  seem  to  have  potential 
for  meeting  this  specification  in  the  next  ten  years. 
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TABLE  6 

1985  PROJECTED  CHARACTERISTICS  1012  BIT  MASS  MEMORIES 


Semiconductor 

Optical 

Electron  Beam 

Magnetic  Bubble 

Type 

CCD 

Bit  by  Bit 

Ion  Implant  Write 

Bubble  Lattice  File 

Volatile/Non- Volatile  Volatile 

Non-Volatile 

Non-Volatile 

Non-Volatile 

Category 

Read/Write 

Archival 

Archival 

Road/Write 

Cell  Sire 

4x4  Micron 

1  Micron 

0.1  Microns 

3x3  Micron 

Chip  Siio 

5  in  Dia.  Wafer 

4  in  x  4  In 

3  in  Dlam.  Wafers  2  In  x  2  In  Wafer 

Chip  Capacity 

2  x  10*  Bits 

10»  Bits 

10"  Bits 

1.6  x  10*  Bits 

Modulo  Capacity 

2  x  10*  Bits 

10»  Bits 

10"  Bits 

3.2  x  10*  Bits 

Number  Modules 

5.000 

100 

10 

3125 

Access  Time 

0.1  to  20  mSoc 

1  Sec 

30  uSec 

80  uSec 

Intrinsic  Data  Bate 

10  Mbps 

10  Mbps 

10  Mbps 

0.78  Mbps 

System  Volume 

14  cu  ft 

27  cu  ft 

3.6  cu  ft 

1.4  cu  ft 

System  Weight 

1570  lbs 

300  lbs 

700  lbs 

157  lbs 

System  Power 

Standby  1  to  10 

KW  2500  Watts 

256  Watts 

102  Watts 

System  Cost/Bit 

0.1  to  1  m  4/Bit 

0.06  m  4/Blt  0.02  m*/Blt 

0.25  m  4/91 1 

System  Cost(Total) 

#1.0  to  10  Million 

$0.6  Million 

$0.2  Million 

$  2.5  Million 

biased  on  the  projections  (3)  above,  we  have  shown  that  physically 
small  solid  state  mass  memories  with  at  least  10  bit  capacity,  access 
time  less  than  100  wSec,  weight  several  hundred  pounds  or  less,  and 
power  no  greater  than  a  few  hundred  watts,  are  feasible  in  the  mid-1980's 
time  period  with  several  emerging  technologies.  These  mass  memories  will 
be  key  components  in  real  time  airborne  data  processing  and  change 
detection  systems.  Integrated  battlefield  surveillance  systems  with 
real  time  data  processing  and  change  detection  may  well  be  combined  with 
low  bandwidth,  secure  communication  links  for  direct,  time  shared  command/ 
control  functions  with  multiple,  low-cost  strike  vehicles  to  counter 
masses  of  heavy,  mobile  armor  in  future  tactical  warfare. 
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SECTION  III 

MEASUREMENTS  ON  MAGNETIC  BUBBLE  MATERIALS 
1.  INTRODUCTION 

The  apparent  direction  for  future  work  on  magnetic  bubble  memories 
is  toward  smaller  magnetic  bubbles.  Magnetic  bubble  materials  to  support 
this  work  must  be  characterized  for  pertinent  bubble  properties  as  rapidly 
as  possible.  Current  practice  is  to  measure  the  magnetic  field  for  bubble 
collapse  (H  coll)  and  the  zero-field  stripewidth  (SW)  in  an  optical  micro¬ 
scope,  to  measure  optical  thickness  (d  opt  =  n  t,  where  n  is  the  optical 
index  of  refraction  and  t  is  the  physical  thickness),  and  to  measure 
magnetic  anisotropy  in-plane/out-of-plane  either  in  an  M-H  looper,  in  a 
microwave  resonance  apparatus,  or  in  a  special  microscope  capable  of  large 
in-plane  magnetic  fields^.  This  measurement  is  limited  to  bubbles  that 
can  be  resolved  in  an  optical  microscope,  that  is,  at  least  three  micro¬ 
meters  in  diameter.  The  thickness  measurement  must  use  an  assumed  (or 
measured,  but  the  measurement  is  seldom  made)  optical  index  of  refraction 
and  requires  a  second  piece  of  apparatus,  an  optical  absorption  spectro¬ 
meter.  The  anisotropy  measurement  requires  an  in-plane/out-of-plane 
magnetic  measurement  and  a  third  piece  of  apparatus.  This  set  of  measure¬ 
ments  fully  characterizes  the  static  magnetic  properties  of  bubble  materials. 

An  improved  technique  for  characterizing  magnetic  bubble  materials 
has  been  described  which  uses  diffraction  from  the  random  array  of  stripes 
at  zero  applied  field  and  at  the  field  for  magnetization  equals  half  the 
saturation  value^.  This  "Laser  Spatial  Filtering"  (LSF)  technique 
eliminates  the  need  for  resolution  of  stripe  edges  and  in  principle, 
reduces  the  optical  characterization  to  a  diffration  measurement,  plus 
the  thickness  measurement  and  anisotropy  measurement.  Unfortunately, 
our  experience  is  that  as  the  bubble  film  gets  thinner  and  thinner  (for 
smaller  and  smaller  bubbles),  the  diffracted  intensity  gets  smaller  and 
smaller  and  becomes  unusable  for  films  thinner  than  about  three  micro¬ 
meters.  Note  that  LSF  characterization  still  requires  three  distinct 
pieces  o'  apparatus  for  static  characterization. 
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In  an  effort  to  simplify  the  measurement  of  static  bubble  materials 

properties,  a  new  technique  was  investigated  for  microwave  measurement  of 

1 9 

the  magnetic  fields  for  bubble  collapse  and  stripout  w.  The  bubble 
material  would  already  be  mounted  for  straightforward  in-plane/out-of-plane 
magnetic  anisotropy  measurements  in  the  same  apparatus.  It  was  felt  that 
the  LSF  technique  would  be  applicable  to  stripewidth  and  saturation 
magnetization  measurements.  The  LSF  was  thought  to  be  capable  of  defect 
mapping  at  that  time,  all  in  the  same  apparatus.  Thickness  was  thought 
to  be  a  quantity  that  could  be  derived  magnetically  from  these  measurements. 
Thus,  complete  static  characterization  could  have  been  done  with  two  pieces 
of  apparatus,  freeing  the  microscope  for  other  work.  Unfortunately,  the 
LSF  could  not  be  made  to  determine  saturation  magnetization  for  films 
thinner  than  three  micrometers.  Additionally  it  has  been  learned  that 
the  LSF  does  not  map  defects  other  than  physical  scratches  and  dust  on 
bubble  films.  This  approach  was  abandoned  as  having  too  many  elements 
not  yielding  useful  data.  Should  an  LSF  technique  for  measuring  saturation 
magnetization  of  small  bubble  films  be  found,  or  should  an  LSF  technique 
for  mapping  only  magnetic  defects  of  small  bubble  films  be  found,  the 
microwave  -  LSF  combination  would  be  worth  reconsidering. 

The  best  combination  of  techniques  for  measurement  of  static  bubble 
properties  was  found  to  be  zero-magnetization  susceptibility  (slope  of  the 
M-H  curve  at  zero  magnetization)  times  thickness,  field  for  bubble  stripout, 
and  field  for  bubble  collapse^’^.  All  these  quantities  can  be  measured 
in  a  suitable  M-H  looper.  The  anisotropy  measurement  can  be  made  in  a 
microwave  resonance  apparatus.  From  these  quantities,  all  of  the  static 
bubble  parameters  can  be  derived.  The  major  portion  of  this  section  will 
be  devoted  to  this  susceptibility/microwave  combination  technique.  Further 
work  is  being  done  to  investigate  the  feasibility  of  using  in-plane/out-of¬ 
plane  susceptibility  to  determine  the  anisotropy,  to  investigate  the  M-H 
looper  "stickiness"  (dependence  on  modulation  field  amplitude)  for  determining 
coercive  field  H^,  and  to  investigate  the  M-H  looper  frequency  dependence 
for  determining  domain  wall  mobility.  In  other  words  two  pieces  of  apparatus 
are  needed  for  static  characterization  now,  and  the  directions  are  evident 
for  reducing  this  to  a  one  piece  apparatus  which  also  perform  the  dynamic 
materials  characterization. 
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2.  MICROWAVE  TECHNIQUE  FOR  CHARACTERIZING  BUBBLE  MATERIALS 

During  a  systematic  investigation  of  the  ferromagnetic  resonance 

absorption  characteristics  of  magnetic  bubble  materials,  a  low- field 

change  in  microwave  susceptibility  was  observed,  as  shown  in  Figure  2, 

in  addition  to  the  g=2  resonance.  This  susceptibility  shift  was  found 

to  be  repeatable  but  was  dependent  on  whether  the  field  was  decreasing 

or  increasing.  The  shifts  were  found  to  be  correlated  with  the  bubble 

collapse  field  or  stripe  collapse  field  for  increasing  field  sweep  and 

1 2 

the  bubble  strip-out  field  for  decreasing  field  sweep 

At  X-band  frequency,  a  change  in  microwave  susceptibility  was  observed 
at  low  fields  for  the  magnetic  field  in  the  out-of-plane  orientation. 

Operating  the  spectrometer  in  the  derivative  presentation  mode,  this 
change  appeared  as  a  peak  with  a  distinctive  nonresonance  lineshape  as 
shown  in  Figures  2-4,  depending  on  the  magnetic  history  of  the  bubble 
material  sample.  If  the  sample  contained  stripes,  we  obtained  the  spectrum 
shown  in  Figure  2  for  either  direction  of  field  sweep  (positive  or  negative 
field).  If  a  large  (5  KOe)  in-plane  magnetic  field  had  been  applied  to  the 
sample,  then  shut  off,  the  metastable  demagnetized  state  was  a  raft  of 
magnetic  bubbles  throughout  the  sample.  If  the  spectrum  was  then  observed 
with  increasing  field  sweep  followed  by  decreasing  field  sweep,  the  spectra 
of  Figures  3  and  4  were  obtained.  Figure  3  represents  bubble  collapse  (up-field 
sweep)  followed  by  "stripe  stripout"  (down  -  field  sweep)  while  Figure  4  presents 
the  case  where  bubble  size  increases  with  increasing  field  with  "negative 
bubble"  collapse  (up-field  sweep)  followed  by  "stripe  stripout"  (down-field 
sweep).  We  observe  that  the  shape  and  position  of  the  susceptibility  change 
with  decreasing  field  sweep  is  indistinguishable  for  the  three  cases  and 
also  if  the  field  was  increased  (in  the  direction  such  that  bubbles  decreased 
in  size  with  increasing  field)  to  within  the  bubble  stability  range,  then 
decreased  past  the  bubble  stripout  field.  Put  another  way,  we  observe 
that  bubble  stripout,  stripe  "stripout",  and  "negative  bubble"  stripout 
are  indistinguishable  in  position  and  lineshape  when  measured  by  low-field 
microwave  susceptibility  changes.  On  the  other  hand,  bubble  collapse, 
stripe  collapse,  and  "negative  bubble"  collapse  have  distinctive  lineshapes 
and  occur  at  different  applied  magnetic  fields.  We  observed  similar 
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Figure  2.  Low-Field  Microwave  Spectrum  of  Stripe  Collapse  and  Stripout 
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Figure  3.  Low-Field  Microwave  Spectrum  of  Bubble  Raft  Collapse  and  Stripout 
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Figure  4.  Low-Field  Microwave  Spectrum  of  "Negative  Bubble"  Collapse  and 
Stri pout 
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susceptibility  changes  at  10,  15,  20  MHz  but  did  not  observe  any  low-field 
susceptibility  changes  at  24  GHz.  The  mechanism  responsible  for  the 
susceptibility  change  at  microwave  frequencies  is,  as  yet,  incompletely 
understood  but  is  felt  to  be  due  to  a  spinwave  scattering  from  the  magnetic 
domain  walls. 

With  the  non-optical  technique  for  measuring  the  bubble  collapse  and 

stripout  field,  it  was  decided  to  implement  a  laser  spatial  filter  to 

11  15 

measure  stnpewidth  and  saturation  magnetization  ’  .  These  measurements, 

when  combined  with  optical  spectrophotometer  measurements  of  film  thickness 
and  microwave  ferromagnetic  resonance  anisotropy  measurements  provide  a 
complete  static  characterization  of  the  bubble  material.  Note  that  in  this 
measurement  sequence,  low- field  microwave  susceptibility  measurement  of 
bubble  collapse  or  stripout  field  is  non-optical,  microwave  ferromagnetic 
resonance  anisotropy  measurements  are  non-optical,  laser  spatical  filter 
measurements  are  optical  but  are  based  on  optical  diffraction,  and  optical 
spectrometer  thickness  measurements  are  based  on  interference  effects. 

In  addition,  the  laser  spatial  filter  was  thought  to  be  capable  of  whole- 
wafer  defect  mapping  which  is  a  measure  of  quality  in  the  growth  procedure. 
Our  experience  with  laser  spatial  filtering  is  reported  in  detail  elsewhere, 
but  will  be  summarized  here. 

Laser  spatial  filtering  is  a  technique  for  diffracting  an  incident 
polarized  laser  beam  from  an  array  of  stripe  domains  in  the  bubble  material 
sample.  When  the  sample  is  preconditioned  to  be  in  a  relaxed  state  at  zero 
applied  magnetic  field,  a  first  order  diffraction  pattern  is  observed  which 
can  be  used  to  determine  the  zero-magnetization  stripe-width.  From  this 
stripewidth,  the  position  of  the  second  order  diffraction  can  be  calculated. 
If  the  second  order  pattern  is  observed,  it  reaches  a  maximum  in  intensity 
when  the  sample  is  magnetized  to  half  of  the  saturation  value.  The  magnetic 
field  to  achieve  this  half-saturation  value  is  noted  and  can  be  used  to 
determine  the  saturation  magnetization.  When  we  sought  to  implement  the 
laser  spatial  filtering  measurement,  we  found  severe  difficulty  with  pre¬ 
conditioning  the  sample  to  achieve  a  reproducible  and  believable  zero- 
magnetization  state  with  four  or  five  micrometer  thickness  bubble  material. 
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This  was  partially  solved  by  applying  a  combined  ac  and  dc  field  with  dc 
value  exceeding  the  collapse  field,  then  decreasing  the  dc  field  to  zero. 

This  procedure  allowed  us  to  determine  the  zero-field  stripewidth  with 
about  2%  reproducibility.  The  spatial  filter  was  then  moved  to  a  calculated 
location  and  the  second-order  diffraction  pattern  focused  into  the  silicon 
PIN  detector. 

The  dc  magnetic  field  on  the  sample  was  increased  until  the  intensity 

was  a  maximum.  For  four  or  five  micrometer  thick  materials,  this  procedure 

was  reproducible  to  within  10%.  The  materials  which  are  of  real  interest 

to  this  effort  are  one  or  two  micrometers  thick  and  the  next  step  was  to 

attempt  to  measure  them.  We  found  that  the  first  order  diffraction  intensity 

was  reduced  by  a  factor  of  about  ten,  which  decreased  the  reproducibility 

in  the  measurement  of  zero-field  stripewidth  to  about  10%.  Repeated 

attempts  to  observe  the  second-order  diffraction  were  completely  unsuccessful 

due  to  extremely  low  intensity.  It  was  finally  concluded  that  laser 

spatial  filtering,  in  its  present  state,  cannot  be  used  for  characterizing 

bubble  materials  thinner  than  about  three  micrometers,  primarily  due  to  low 

intensity.  At  a  late  stage  in  this  investigation  it  was  learned  that  the 

laser  spatial  filter,  in  its  present  state,  cannot  be  used  to  map  magnetic 

defects  because  physical  scratches  and  dust  on  the  surface  completely 

1  5 

obscure  any  magnetic  defects  that  might  be  present. 

The  microwave/LSF  approach  to  characterizing  one  and  two  micrometer 
thick  magnetic  bubble  materials  was  finally  abandoned,  not  because  the 
microwave  measurement  failed,  but  because  the  laser  spatial  filtering 
technique  had  too  many  components  which  were  not  yielding  useful  data  on 
small  bubble  materials.  The  original  hope,  that  zero-field  stripewidth 
and  saturation  magnetization  could  be  used  to  derive  an  effective  magnetic 
thickness,  was  not  carried  further  due  to  the  apparent  impossibility  of 
obtaining  reliable  experimental  data.  Before  the  approach  was  abandoned, 
a  computer  program  was  written  to  derive  the  static  magnetic  bubble 
parameters  from  bubble  collapse  field,  zero-field  stripewidth,  and  physical 
thickness  using  the  theory  of  Thiele. ^  This  computer  program  is  reproduced 
in  the  Appendix  with  instructions  for  use  and  sample  calculations.  This 
program  can  also  derive  the  static  magnetic  bubble  parameters  from  laser 


19 


A  FAL-TR- 79-11 28 


spatial  filter  data:  field  for  magnetization  equal  to  half  the  saturation 

value,  zero-field  stripewidth,  and  physical  thickness.  In  this  computer 

program,  the  exact  theory  of  Thiele^  was  used,  rather  than  the  commonly 

applied  approximation  of  Callen  and  Joseph,^  to  explore  the  accuracy  of 

the  approximate  calculation  for  small  bubble  materials.  Cerlcinly  the 

exact  theory  program  runs  very  rapidly  and  does  not  require  much  computer 

memory;  the  advantages  of  the  approximate  theory  seem  very  small.  In 

fact,  even  the  exact  theory  (which  does  not  include  the  effects  of 

anisotropy)  seems  to  introduce  errors  into  the  calculation  of  materials 
1 5 

parameters  and  further  approximations  may  worsen  the  errors. 

3.  SUSCEPTIBILITY  TECHNIQUES  FOR  CHARACTERIZING  BUBBLE  MATERIALS 

A  technique  for  bubble  materials  characterization  is  described  which 
uses  low  frequency  magnetic  susceptibility  and  bubble  collapse  and 
stripout  measurements.  This  method  is  presented  as  an  alternative  to  the 
microwave  and  laser  spatial  filtering  technique.  It  can  be  automated  and, 
since  domains  are  not  detected  optically,  there  are  no  resolution 
problems.  The  technique  can  also  be  used  on  opaque  magnetic  metal  films. 

The  low  frequency  susceptibility  characterization  determines  the 
zero-magnetization  stripwidth,  characteric  length,  saturation  magneti¬ 
zation,  and  thickness  from  the  zero-magnetization  susceptibility 

=X„,the  bubble  stripout  field  H  ,  and  the  bubble  collapse  H 
dH  o  r  so  r  coll 

The  thickness  is  a  magnetically  determined  effective  value,  which  could 

differ  from  the  physical  thickness  due  to  dead  layers  or  other  nonmagnetic 

layers  deposited.  Henry  has  pointed  out  that  the  difference  between 

magnetic  thickness  and  optical  thickness  will  become  significant  for 

films  below  one  micrometer  thickness. 
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1  O 

The  theory  of  Kooy  and  Enz  leads  to  the  following  two  equations 
which  interrelate  the  static  properties  of  bubble  domain  materials: 


0  =  +  —x—  I  — s-  sin  {n  tt  (1  +  JJ-)}  {1  -  exp(-  2*nh)} 


4ttM, 


S  tt  h  n=l  n 


0) 


t  _  d  Z  .2  rniT  n  ,  M  rl  ,  2irnh>  <  2irnh,, 

F1J  *  sin  {2“(1+MC)}  {1  -  0  +  —  )exp(-  -*-)> 

it  h  n=l  S 


(2) 


Here  is  the  saturation  magnetization,  M  is  the  average  magnetization 
of  the  film,  d  is  the  domain  period,  h  is  the  film  thickness,  and  H  is  a 
magnetic  field  applied  perpendicular  to  the  film's  surface.  The  derived 
quantity  l  is  the  characteristic  length,  a  materials  constant. 

The  magnetic  susceptibility  jjjj-  =  X  can  be  obtained  by  solving 
equations  (1)  and  (2)  numerically  and  evaluation  the  slopes  of  the 
resulting  M  vs.  H  curves.  However,  our  interest  is  only  in  the  zero- 

p 

magnetization  susceptibility  X  ,  and  the  corresponding  value  of  After 
some  algebraic  manipulations,  the  following  equations  are  obtained  for 
M  =  0: 


and 


4"*0  ■  <>  +iff 


n=l 


[1  -  exp  (  - 


l  _ 
h 


4w 


TT  h 


n(odd)  n" 


K  [i 


0  *  «P  (-  “?)] 


where  w  is  the  zero-magnetization  stripewidth. 


(3) 


(4) 
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The  calculated  values  SW/h,  4wXq,  ^£°  ,  and  H$0/Hcoll  are 

n  S  S 

plotted  as  functions  of  jj-  in  Figures  5  and  6.  A  computerized  analysis 
procedure  using  tables  of  values  of  the  quantities  in  Figures  5  and  6 
has  been  written  to  analyze  experimental  data  and  calculate  the  magnetic 
bubble  materials  parameters.  This  program  is  listed  in  the  Appendix  with 
examples  of  its  use.  At  least  three  out-out-plane  experimental  quantities 
must  be  measured,  including  one  each  magnetic  and  linear  quantities.  Here 
zero-magnetization  susceptibility  4ttX  is  considered  a  magnetic  quantity 
(but  an  additional  magnetic  quantity  must  be  measured  to  completely 
characterize  the  material),  while  susceptibility  times  thickness  4-rrXgh 
(the  actual  experimental  quantity)  is  considered  a  linear  parameter.  This 
program  calculates  a  table  of  bubble  diameters  vs  applied  field  if  charac¬ 
terization  is  complete,  and  will  also  perform  theoretical  calculations  of 
all  the  other  parameters  if  a  complete  set  of  materials  parameters  is 
input.  This  program  is  listed  in  the  Appendix  together  with  instructions 
for  use  and  sample  calculations. 
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Figure  5.  Plot  of  Rationalized  Zero-Magnetization  Susceptibility  (4tXq) 

and  Zero-Magnetization  Stripewidth  by  Thickness  (SW/h)  Against 
Characteristic  Length  by  Thickness  (£/h) 


!  /h 


Figure  6.  Plot  of  Bubble  Collapse  Field  by  Saturation  Magnetization, 

Bubble  Stripout  Field  by  Saturation  Magnetization,  and  Bubble 
Stripout  Field  by  Bubble  Collapse  Field  Against  Characteristic 
Length  by  Thickness 
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SECTION  IV 

CONCLUSIONS  AND  RECOMMENDATIONS 

Program  BUBCHAR  is  written  in  CDC  FORTRAN  EXTENDED,  but  an  effort 
has  been  made  to  use  FORTRAN  IV  coding  wherever  possible.  Thus  the 
program  should  run  on  most  computers  with  few  changes. 

There  appears  to  be  small  advantage  to  using  the  Callen  and  Josephs 

approximation^ 7  in  bubble  materials  calculations.  The  ratios  in  the  data 

16  18 
statements  were  computed  from  Thiele  theory  or  Kooy  and  Enz  theory  . 

In  fact,  two  instructions  are  included  to  unraval  the  saturation  magneti¬ 
zation  calculated  using  the  Callen  and  Josephs  approximation,  then  re¬ 
calculate  the  materials  parameters  according  to  Thiele  theory^.  The 
major  inaccuracy , in  analysis  is  the  ad  hoc  introduction  of  anisotropy, 
after  calculating  all  other  parameters  assuming  that  the  material  is 
isotropic. 

The  computer  program  BUBCHAR  is  listed  in  the  Appendix  following, 
and  if  the  author  is  sent  a  blank  half-inch  magnetic  computer  tape,  he 
will  write  card  images  of  the  program  on  the  tape  and  return  it.  Be  sure 
to  specify  density,  number  of  tracks,  and  any  standards  or  conventions 
that  would  make  the  tape  easier  to  use.  As  a  poor  alternative,  the  author 
can  provide  a  card  copy  of  the  program  (Be  sure  to  specify  the  keypunch 
standard  needed.  Only  "029"  keypunch  standard  cards  can  be  interpreted 
at  this  facility). 
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APPENDIX 


A.  INTRODUCTION 

A  CDC  Fortran  Extended  computer  program  has  been  written  primarily 
to  analyze  sets  of  experimental  data  and  calculate  static  bubble  materials 
parameters.  Three  out-of-plane  measurements  are  needed,  plus  one  in-plane 
measurement  to  specify  the  anisotropy.  The  computer  then  calculates 
interpolations  and  ratios  analogously  to  the  graphical  procedure  described 
by  Thiele^  to  determine  the  bubble  materials  static  parameters.  The 
program  will  also  accept  a  complete  set  of  bubble  materials  parameters 
and  calculate  values  of  all  the  experimental  and  materials  parameters. 

If  a  sufficient  set  of  out-of-plane  parameters  is  input,  the  program 
calculates  a  set  of  stable  bubble  diameters  and  applied  field  values. 

If  no  in-plane  information  is  input,  the  program  simply  does  not  calculate 
anisotropy  field,  anisotropy  energy  density  or  q-factor. 

B.  INSTRUCTIONS  FOR  INPUT 

Program  BUBCHAR  is  written  to  accept  data  cards  consisting  of  an 

instruction  (columns  1-10)  followed  by  a  number  (containing  a  decimal 

point)  in  columns  11-20.  Non-instruction  comment  cards  may  be  used 

without  restrictions.  The  permissible  instructions  are  listed  and 

defined  in  Table  A1 .  As  soon  as  three  out-of-plane  parameters  are 

input  (including  at  least  one  each  magnetic  and  linear  quantities),  the 

program  calculates  all  the  other  quantities  from  Thiele  Theory^  or  Kooy 
1 8 

and  Enz  Theory  .  An  in-plane  parameter  may  be  input  anywhere  in  the 
data  deck  or  may  be  omitted. 
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TABLE  A-l 


INSTRUCTIONS  FOR  INPUT  DATA:  PROGRAM  BUBCHAR 


Instruction 

Classification 

Definition 

HC 

out  of  plane  magnetic 

bubble  collapse  field  H  (Oe) 

HCOLL 

•1 

II  C 

hcollapse 

•I 

II 

HSO 

ft 

bubble  stripe  out  field  H  (Oe) 

li  so 

HSC 

•1 

HSTRIPEOUT 

II 

II 

T 

out  of  plane  linear 

bubble  material  thickness  h  (microns) 

THICK 

II 

II 

THICKNESS 

•1 

It 

•1 

II 

H 

X0H 

II 

aero-magnetization  susceptibility  times  thickness 

CHI0H 

M 

II 

SW 

•1 

II 

zero-magnetization  stripewidth  (microns) 

It 

s 

STRIPEWIDT 

•I 

II 

STRIPWIDTH 

II 

II 

MAG 

out  of  plane  magnetic 

saturation  magnetization  (gauss) 

MS 

II 

II 

4. PI. MS 

•1 

•1 

M 

II 

II 

L 

out  of  plane  linear 

characteristic  length  P  (microns) 

CHARLEN^TH 

•1 

II 

X0 

Our  of  plane  magnetic 

zero-magnetization  susceptibility  (Gauss/Oe) 

CH10 

•t 

II 

Q 

in  plane  magnetic 

anisotropy  quality  factor 

HU 

•1 

anisotropy  field  (Oe) 

HA 

«• 

ii 

HK 

•I 

ii 

KU 

•1 

anisotropy  energy  density  (ergs/cm^) 

4.PI.M 

out  of  plane  magnetic 

saturation  magnetization  calculated  on  Callen &  Joseph  appro* 

CJM 

•• 

II 

H(l/2) 

•i 

field  for  half-saturation  magnetization  (Oe) 

Hl/2 

•• 

it 

H|X 

ii 

ii 

HHALF 

•i 

ii 

LABEL 

Alphanumeric  identifier 

(up  to  70  characters) 

TITLE 

Any  Other 

M 

Comments  (printed  as  they  are  punched)  Do  not  use  just  before  an  in-plane  card. 
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C.  EXAMPLES  OF  USE  OF  BUBCHAR 

1.  Example  1: 

Input:  H 
SW 

HCOLL 

HK 


Output:  Fig  A-l 

2.  Example  2: 

Input:  H 

CHI0H 

HSO 

HK 

Output:  Fig  A-2 

3.  Example  3: 

Input:  H 

4. PI. MS 

L 

Q 

Output:  Fig  A- 3 

4.  Example  4: 

Input :  H 
X0 
SW 

Output:  Fig  A-4 


D.  Listing  of  Program  BUBCHAR 
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PROGRAM  BUBCHAR 


H  2.40000 
SW  1.800CO 
HCOLL  409.00000 


H  =  . 24C  0  0  3E*3 1 ,  L  =  .  169471E*-00,  OHIO  =  .1478616+01,  SW  =  .183D0OE+01 

4. PI. MS  =  .673001E+03,  HS  =  .3301226+03,  HC  =  .  4  C  90  0  ClE  +  0  3,  BULK  WALL  ENERGY  =  .61082SE+C0 

Hl/2  (FIELO  FOR  HALF-SATURATION  MAGNETIZATION)  =  .219570E+03 


BUBBLE  DIAMETER  {MICRONS) 

APPLIED  FIELO  (0E) 

1.12 

406.52 

1.29 

400.56 

1.46 

392,75 

1.63 

384.06 

1.  33 

374.93 

1.97 

365.66 

2.14 

356.44 

2.33 

347.41 

2.47 

338.62 

BUBBLE  COLLAPSE  DIAMETER  =  .95  MICRONS 

BUBBLE  COLLAPSE  FIELO  -  409.00  OERSTEDS 

BUBBLE  RUNOUT  DIAMETER  =  2.64  MICRONS 

BUBBLE  PUNOUT  FIELD  =  330.12  OERSTEOS 


HK  4172.00000 


MAXIMUM  TWIST  ANGLE . .  88.4342  DEGREES  (  1.54  RAD.) 

OOMAIN  WALL  ENERGY . .  .5858S5E+00  ERGS/SQ.CM. 

WALL  THICKNESS . .  .-.48362E-01  MICRONS 


Q  =  6.199,  KU  =  .1U717E  +  06,  HK  =  .417200E+C4,  AND  EXCHANGE  CONSTANT  =  .2C8735E-06 


Figure  A-l  Output  from  computer  program  BUBCHAR  for  Example  1.  This  calculation 
is  for  the  standard  characterization  measurements  with  thickness  H  = 
2.40  micrometers,  zero-magnetization  stripewidth  SW  =  1.80  micrometers, 
bubble  collapse  field  HCOLL  =  409.0  Oersteds,  and  anisotropy  field  HK  = 
4172.0  Oersteds. 
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PPOGRAM  8U8CHAR 


H  2.40000 
CHICH  3.30030 
HSO  299.60000 


H  =  .2L0003E+31,  L  =  .122C33E  +  00,  OHIO  =  .  1375 C C E  +  G  1 ,  SW  =  .1497722+01 

4. PI. MS  =  .5453552+03,  HS  =  .29960CE+C3,  HC  =  .361636E+QJ,  BULK  MALL  ENERGY  =  .289567E+00 

K 1/2  (FIELD  FOR  HALF-SATURATION  MAGNETIZATION)  =  .1931282+03 


BUBBLE  DIAMETER  (MICRONS) 

APPLIED  FIELD  (0E> 

.92 

359.79 

1.35 

355.34 

1.20 

349.42 

1.34 

3**2. 65 

1.43 

335.48 

1.61 

323.18 

1.75 

320.87 

1.39 

313.63 

2.03 

306.53 

BUBBLE  COLLAPSE  DIAMETER 
BUBBLE  COLLAPSE  FIELD 
BUBBLE  RUNOUT  DIAMETER 
BUBBLE  RUNOUT  FIELD 


.78  MICRONS 
361.64  OERSTEDS 
2.17  MICRONS 
299.60  OERSTEOS 


HK 


1273. 00000 


MAXIMUM  THIST  ANGLE . =  87.3025  DEGREES  (  1.52  RAO.) 

DOMAIN  HALL  ENERGY . =  .262602E+00  ERGS/SQ.CM. 

HALL  THICKNESS . =  . 90 364SE-0 1  MICRONS 

Q  =  2.  332,  KU  =  .2765 32E  +  0 5 ,  HK  =  .  127300E  +  C4,  ANO  EXCHANGE  CONSTANT  =  .  189510E-06 


Figure  A-2  Output  from  computer  program  BUBCHAR  for  Example  2.  This  calculation 
is  for  the  suggested  optimum  characterization  measurements  with  the 
thickness  H  =  2.40  micrometers,  zero-magnetization  susceptibility 
times  thickness  CHIOH  =  3.30  EMU-micrometers ,  bubble  stripout  field 
HSO  =  299.6  Oersteds,  and  anisotropy  field  HK  =  1273.  Oersteds. 
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PROGRAM  BUBCHAR 


H  3.10000 
4. PI. MS  537.50000 
L  .09740 


H  =  .310000E*01f  L  =  .97400CE-01,  OHIO  =  .127063E4C1,  SH  =  .149998Ef01 

4. PI. MS  =  .537500E*03,  HS  =  . 336132E+03,  HC  =  .392624E+03,  BULK  WALL  ENERGY  =  .223927E+00 

Hl/2  (FIELO  FOR  HALF-SATURATION  MAGNETIZATION)  =  .207544E4G3 


BUBBLE  DIAMETER  (MICRONS) 

APPLIED  FIELO  (OE) 

.91 

390.91 

1.04 

387.01 

1.18 

381.83 

1.32 

375. 89 

1.45 

369.49 

1.59 

362.84 

1.73 

356.09 

1. 86 

349.34 

2.03 

342.69 

BUBBLE  COLLAPSE  DIAMETER  =  .77  MICRONS 

3UBBLE  COLLAPSE  FIELO  =  392.64  OERSTEDS 

CUBDLE  RUNOUT  DIAMETER  =  2.14  MICRONS 

BUBBLE  RUNOUT  FIELO  =  336.15  OERSTEDS 


Q  3.20000 


MAXIMUM  TWIST  ANGLE . =  83.522C  OEGREES  <  1.55  RAO.) 

DOMAIN  WALL  ENERGY . =  .214858E+03  ERGS/SQ.CM. 

WALL  THICKNESS . =  •  498552E- 0 1  MICRONS 


Q  =  3.  200  ,  KU  =  .  3678  47E*  0 5 ,  MK  =  .  1720  30E  +  04 ,  AND  EXCHANGE  CONSTANT  =  .851971E-07 


Figure  A-3  Output  from  computer  program  BUBCHAR  for  Example  3.  This  calculation 
is  a  theoretical  calculation  for  a  bubble  material  with  thickness  H  = 
3.10  micrometers,  rationalized  saturation  magnetization  4. PI .MS  =  537.5 
Gauss,  cnaracteristi^  length  L  =  0.0974  micrometers,  and  anisotropy 
quality  factor  Q  =  3.2. 
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PROGRAM  3U8CHAR 

H  1.50030 

XO  1.40930 

SM  1.03330 

XC  WAS  THE  ONLY  MAGNETIC  PARAMETER  INPUT 

XC  =  .  14C937E*- jl,  H  =  .  15  3 0 j OE  +  ul »  SW  =  .103Q0CE»C1,  AND  L  =  .860438E-01. 

ALL  That  can  be  said  ABOUT  the  OTHER  magnetic  parameters  is 
HC/“.  PI .  MS  =  .o431 15E  ♦  0  AND  HS/4.PI.MS  =  .827942E*00. 


Figure  A-4  Output  from  computer  program  BUBCHAR  for  Example  4.  Note  that  the 

characterization  is  incomplete,  as  zero-maanetization  susceptibility 
(or  zero-magnetization  susceptibility  times  thickness)  is  insufficient 
to  specify  the  magnetic  parameters  of  a  material.  Input  was  thickness 
H  =  1.50  micrometers,  zero-magnetization  susceptibility  XO  =  1.409  EMU, 
and  zero-magnetization  stripewidth  SW  =  1.00  micrometers. 
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PROGRAM  BU8CHAR(I NPUT=1018, OUT PUT=1C1B,TAPE5=IN PUT, TAPE 6= OUTPUT)  8U0CH  1 

01  HENS  I  ON  PHI 1(123),  EOASH ( 1 20 ) , RAT  1 0  ( 1 20 >  BU8CH  2 

DIMENSION  HCBYMS(4G> ,HC8YHS(40)  BUBCH  3 

DIMENSION  HSBYHSI4C)  , WDV H  (4 0 )  , CHI  0 ( 4  C  >  , LABEL ( 8)  BUBCH  4 

DIMENSION  INSTR<35> ,CARQ (8) ,KMAG<2) ,KLEN<2) ,01 (56)  BUBCH  5 

DIMENSION  HHALFUO)  BUBCH  6 

REAL  L,MS,L0H,L0VH(4C) ,KU  BUBCH  7 

COMMON  /AA/  Q  BUBCH  8 

DATA  DI/56*1./,LABEL/8*1H  /  BUBCH  9 

BUBCH  ID 

1=  1  2  345  67  3U8CH  11 

BUBCH  12 

DATA  INSTR/ '*HC"» "HCOLL",  "HCOLL APSE", “HSO", "HSC " , "HSTR I POU T" , "T",  BUBCH  13 

BUBCH  14 

8  9  10  11  12  13  14  15  BUBCH  15 

BUBCH  16 

1  "THICK", "THICKNESS", "H" , "X  OH" ,"CHI 0 H", "SW" , "S" ,"STRI PE  HI DT" ,  BUBCH  17 

BUBCH  18 

16  17  18  19  2  G  21  22  23  BUBCH  19 

BUBCH  20 

2  "STRIPWIDTH" ,  "M  AG",  "MS"  ,  "4  .  PI .  MS"  ,  “M“  ,  "L",  "CHAR1.  ENGT  H“  ,"XG",  BUBCH  21 

BUBCH  22 

24  25  '26  27  28  29  30  31  32  33  BUBCH  23 

BUBCH  24 

3  "CHI C " , "Q" > "HU" , "HA", "HK", "KU", "4. PI . M" ,"C JM" , "H ( 1/2 ) " ,"H1 /2" ,  BUBCH  25 

C  BUBCH  26 

C  34  35  BUBCH  27 

4  "H 12" , "HHAL F“/  BUBCH  28 

C  BUBCH  29 

0DATA  PHU/6*0.  ,.  01, .03, .  05,.  07,. 1,  .4,  .58, .9,1.  14,1.  26,1.33,1.  375,  BUBCH  30 

1  1.  4,1.42,1.435,1.48,1.463,1.5,1.52,1.525,1.531,1.537,1.54,1.545,  BUBCH  31 

2  6*0. ,.  01, .  03, .05, .  07, .  08, .  35, .5, .  865, 1.08, 1.23, 1.31,  BUBCH  32 

3  1. 35, 1.38, 1.41, 1.42, 1.47, 1.48 5,1.?, 1.52, 1.525,1. 531, 1.  537, 1.  54,  BUBCH  33 

4  1. 545,6*0. , . 01, . 33,. 04, . 065 , . 07, . 27 , .4,. 76, 1. , 1. 14,1. 26, 1. 32 ,  BUBCH  34 

5  1. 35, 1 .4,1 .41,1. 48,1.485 , 1. 5,1. 52,1. 525,1. 531 , 1.  537, 1. 6-t  BUBCH  35 

6  ,1.545,6*0. ,. 00 1,. 009,. 01, . 02, . 03, . 13,. 2,. 42,  65 , . 84 , . 99 , 1 . 1 ,  BUBCH  36 

7  1.2,1.25,1.3,1.45,1.47,1.5,1.52,1.525,1.531,1  537,1.54,1.545/  BUBCH  37 
0DATA  EDASH/. 825, . 94, .98,1. ,1. 2,1.  24,1.28,1.  32, 1. 36,1.  4, 1. 45,1.5,  BUBCH  38 

1  1.5,1.45,1.4,1.37,1.372,1.4,1.418,1.44,1.47,1.59,1.6,1.43,1.52,  BUBCH  39 

2  1.56,1.58,1.8,1.82,1.325,1.08,1.18,1.22,1.28,1.34,1.4,1.43,1.47,  BUBCH  40 


3  1. 5, 1. 52, 1. 55,1. 565, 1. 58 , 1. 52,1.  43, 1.44,1.  43, 1.4  5,1.  46, 1,48, 1.52, BUBCH  41 

4  1.61,1.61,1.49,1.53,1.57,1.59,1.31,1.83,1.835,1.5,1.525,1.55,  BUBCH  42 

5  1.575,1.6,1.625,1.65,1.675,1.7,1.72,1.74,1.75,1.749,1.71,1.7,  BUBCH  43 

6  1.625,1.632,1.65,1.67,1.675,1.69,1.72,1.75,1.79,1.805,1.32,  BUBCH  44 

7  1.  335,1.  842,1.85,1.855,1.63,1.84,1.85,1.855,1  06, 1.87,1  88,1.9,  BUBCH  45 

8  1. 90 3, 1.90 7, 1.91, 1.91 5, 1.92, 1.91, 1.9, 1.895, 1.  3 94, 1.89,1.  3a 5 , 1 . 883DUBCH  46 

9  ,1.83,1.355,1.59,1.895,1.9,1.905,1.91,1-913,1  916,1.919/  BUBCH  47 

OOATA  RATIO/. 35,.  365,.  375  ,.  <*1, . 44,. 465 ,.49, .52, .54, .56, . 538,. 7  ,  BUBCH  48 

1  .’2, .  755, .  74,. 725, .73, .  725, .73, .74, .  75, .795, .  8, .84,.  86,  87,.  83  BUBCH  49 

2  ,.  89, .  895, .9, .43, .49, .5, .527, .55, .565, .59, .618, .63, .645, .66,  3UBCH  5 0 

3  .  75, .  76, .  778, .755, .74, .738, .738, .746, .  75, .76, .  805, .81, .85, .87,  BUBCH  51 

4  .  38, .  39, .9, .905, .91, .  725, .  73, .735, .75, .76, .73, .79, .797,  8, .803,  BUBCH  52 

5  .  815, .  85, .  86, .86, .85, .83, .82, .82, .82, .82, .83, -345, .855, .  88, *  9,  3UBC H  53 
6-  .908, .  92,. 928, .  928, .  929, .903, .  909, .912, .  915, .918, .92, .  923, .93,  BUBCH  54 

7  .) 35, .  9  38,. 94,. 95,.  955, .  96, . 96,. 95, .  945,. 943, .  93,.  93, .  93,. 92 , .  92, BUBCH  55 

8  .92, .94, .945, .95, .955, .  957, .959/  BUBCH  56 

OATA  PI4/12. 5663702/  BUBCH  57 
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OATA  (LOVHIII, 1  =  1,401/0.. .001. .005,  . 3 1 , .  0 2, .  0  3,  .  0  4,  .  G  5, .  0  6,  .  0  7,  3U3CH  53 

1  .08, .09, .1,. 11 , .12, . 13, .14, .15, . 16, .17, . IS , .19, .2, .21, .22, .23,  BUBCH  59 

2  .24, .25, .26, .27, .28, .29, .3, .31, .32, .33, .34, .35, .36, .37/  9U9CH  60 

DATA  (WOVH(I),I  =  l,43»/  .  0  3  0  0  C  0  6  +  0  0  ,  .  3585446-,.  1  ,  .1919646*00.  3U3CH61 

1  •  27 1 49  3E  +  00 ,  .3844005+03,  .4724416*00,  .548706E+0;,  .  61 8242c*  0  3  ,  3  J9CH  62 

2  . 683589E+30 ,  .7462265*00,  . 80 7 0 92=  ♦0 0 ,  .8668436+01,  . 92 594 75* 0 i  . 3UBC H  63 

3  .984762EOO,  .1043576+31,  .1102616+cl,  .1162356*01,  .  1 222  396*  3 1  ,  BUBCH  6  4 

4  .  1282845*01,  .  1 3444 6E *0 1  ,  .  140 7055* 3 1 .  .147.726*31,  .  1 5 35 5 8  -  + ; l , 9 JBC H  65 

5  .16C17?E*J1,  .  1669256*01 ,  .  17 38 25 £♦. 1 ,  .18C381E+31,  . 1881C25+ I  1, 3UBCH  66 

6  .  195  495E*  01,  .2030725*01,  .21u8396*ul,  .  2  1  88  0  65*  J  1 ,  . 2269025* : 1 , BUBCH  67 

7  .235375E*J1,  .2439945*01,  .2528495*01,  .2629485*01,  . 27 23 32c* j 1 , 3U3C H  6 8 

8  .281921E*01,  .2918125*01/  BUBCH  69 

OATA  (CHIOtl) ,1  =  1,4.)/  .  10  00  L CE*01  ,  .1039375*.  1,  .  1092525  +  0  1,  BUBCH  70 

1  .1136075*01,  .1204125*01,  .1262686*01,  .131741c+01,  . 1 3 70 446* 01 , 3U3CH  71 

2  •  142289E+01,  .  1475386*01  ,  .  1528325*0  1  ,  .1582322*0  1  ,  .  1 6 366 3c *  3 1 , B JBC H  7? 

3  .  1692495*01,  .  1749595*01,  .  1828096*0  1,  .1868115*0  1,  .  1929 7 7£* 3 1 ,  BUBCH  73 

4  .1993145*01,  .2058315*01  . 2 12 5 3 8 5*. 1 ,  .2194415*01,  . 22 655 25* 0 1 , 3 J3C H  74 

5  .233876E*01,  .2414235*01,  .2492215*01,  .2572175*01,  . 26 548 2c+ 2 1 , 9U3C H  75 

6  .2740025*01.  .2827885*31  ,  .  291 8 4 7C ♦.  1  ,  .  3011  905*0  1,  .  31  0S25E* 3 1 , 3J3CH  76 

7  .3207656*01,  .  3310155*31,  .  3415886*1.  1,  .3524945  +  3  1,  . 36 374 46* 0 1 , BUBCH  77 

8  .375348E+01,  .3873195*01/  BUBCH  78 

DATA  (HSBYMSCI) , 1=1,40) /  .luOCOOE+Ol,  .9133165*30,  .3275755*00,  BUBCH  79 

1  .767021F  +  00  ,  .68815  95  +  03  ,  .  632214-+.  3  ,  .5883245+0C,  .  55  IS  8  9c+ 0  3 , 3'JBCH  80 

2  .51998 95 +0J,  .4921275*00,  .4672056*00,  .444614F+0C,  . 42 39925+ 0 3 , BUBCH  81 

3  .4049786*00,  .  3873356  +  03  ,  .  37 ; 9 73  ' *0  0  ,  .3556265  +  00  ,  . 34 122 1 £+ 03  , 8UBC H  02 

4  .327 69 85+0 u,  .3149345*00,  .3028725*10,  .2914006*00,  . 28 C 5 395+ 0 0 , B JBC H  83 

5  .2702276*00,  .2633655+00,  . 250 982 r *c 0 ,  .2419745+0C,  . 2334 366+ 03 , 3UBCH  84 

6  .2252305+00,  .217388E+G0,  .2098406+00,  .20261 26+0l,  . 1 9 57 06E* C 3 , BUBCH  8 5 

7  .1890745+00,  .1826836*00,  .1765236+CO,  .1706126*03,  . 1 649 455+ 30 , BU3CH  86 

8  .1594756*30,  .1541955+03/  3U3CH  87 

DATA  (HCBYMS(I)  ,  1  =  1, 40/0.  10CCGo6*3  1,  .93239.6*00  ,  .  8  88  86  95  +  0  ,  BUTCH  88 

1  .84401 6F  +  00,  .  7826106  +  00  ,  .736317-+03,  .698.7  55  +  00  ,  .  6650  645+  3  -3  ,  BUBCH  89 

2  .6357715+03,  .6092726*03,  .5850035*. 3,  .562550E+3C,  . 54 162 3c+ 0 3 , 3 JBCH  93 

3  .5219955*00,  .5035135*00,  .4860555*. 0,  .4695405+00,  . 45 3361 5+ 0 3 , 3U3C H  91 

4  .4389255*00,  .4246695+00,  .4110335+..,  .3979635*03,  . 38 54 1 7c + 3 . . BUBCH  92 

5  .3733575*00,  .3617545+03,  . 350586"+. 0,  .339331E+00,  . 32 94545+ 0 3 , BUBCH  93 

6  .3194315*00,  .3G9743E+G0,  . 300 3825*. 0 ,  .2913456*00,  . 2825136+ 03 ,(1U3CH  94 

7  .274168E+00,  .2660026+00,  . 2530 946*. 3 ,  .2504305*00,  . 24299 96* 0 0 , BUBCH  95 

8  .2357895*00,  .223788E+00/  3UBCH  96 

DATA  (HCBYHS(I) , I = 1 , 40 ) /0 . lu 0 0C 36+ 0 1 ,  .1020875*31,  .1074065*01,  BUBCH  97 

1  .11O  038E*O1,  .  11372 5E +  01,  .1164665  +  0  1  ,  .1186555  +  0  1,  .  120 55 15* 5 1  , O J 3C H  9 8 

2  .1222666*01,  .1238046+01,  . 1252135*. 1,  .1265266*01,  . 1 2 77 446+ 0 1 , 9U 3C H  99 

3  .128895E  +  01,  .  1299945*01,  .  1313226*.  1,  .132.325+01,  . 1 3 3 j 1 lc* 3 1 , 3JBCH1 0 0 

4  .1339425*01,  .1348445*01,  . 1357126*. 1,  .1365696*01,  . 1 3738 45* 01 , 9U3CH1G 1 

5  .138 164E*0 1 ,  .1339415+01.  .1396365+31,  .1404415*01,  . 14 11325+ 0 1 , BJ3CH1C ? 

6  .1418245*01,  .1424845*31,  . 1431485*. 1,  .1437956*01,  . 1444 0 76* j 1 , BUBCH1 0 3 

7  .1450065*01,  .1456085*01,  .1462105*01,  .1467836*01,  . 1 4 732 16* C 1 , 3 JBCH1 0 4 

8  .147853E+01,  .1483765+01/  3UBCH1C5 

OATA  ( HHALF ( I ) , 1=1 ,40 > /  .4936185+00.  ,4798435+00.  .4549306+00,  3UBCH1C6 

1  .4362436*00,  .4C96085+00,  .3888166*j0,  .3709945*00,  . 3550 556* 00 , BUBCH1 C 7 

2  .3405115*00,  .3270506*00,  . 31448 0 5*. 3 ,  .3c266SE*0G,  . 29 15 156* 0 0 , OUBCHl 0 8 

3  .2839-*9E+0O,  .  2709075  +  00,  .2613446*c  0  ,  .2522185  +  0  0,  .  24  34  9  66*  0  0  .  B'JBC  HI  0  9 

4  .2351436*00,  .22715:6+00,  .2194786*00,  .2121135+00,  . 2C 53 376* J 0 , 3U3CH1 1 C 

5  .1982346  +  00,  .  1916396*00,  .  1  853  8  S-‘*u  J  ,  .  1  79320E*0C,  .  1  7  34  7  35+  0  0  ,  BU3CH1 1 1 

6  .1678385*00,  .162403E  +  00,  .  1571606*0  0,  .1521:05*00  ,  .  1 4 72 1 7c+ 0. . BU3C HI  1 2 

7  .  1425026*00,  .  1379475*00,-  .  1 3 354 95  * o 0  ,  .  129298E  +  0C,  . 1 2 51 91 6* 0 3 , BJ3CH1 1 3 


8  .121 22 16+00,  .1173835*00/ 


BU3CH114 
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9  HS=0. 

10  CONTINUE 

IFIMS.EO.C.)  WRITE(6,5000) 

5000  FORMAT  (“1  PROGRAM  BJBCHAR“//> 

NMAG=Q 

NLEN=0 

KIPM=0 

1  READ(5,960>  INS 

950  FORMAT  (A1G ) 

IF(EOF  (51 )  100,15 

15  CONTINUE 
BACKSPACE  5 
DO  2  1=1,35 

IF ( INSTR(I) .  £Q. I  NS )  GO  TO  3 

2  CONTINUE 
READ(5,961)  CARO 

951  FORMAT  (8  A 1  0  ) 

IF (EOF  (5  I )  100,25 
25  CONTINUE 

WRITE  16,  1060)  CARO 
1 C  50  FORMAT  (IX, 8A10) 

IF(INS, EQ." TITLE". OR. INS. E3." LABEL")  GO  TO  85 
GO  TO  1 

3  REA 0(5, 962)  IN5.A1 

952  FORMAT(A10,F1Ci.C) 

I F ( EOF  (5 ) I  1JG.35 

35  CONTINUE 

WRITE ( 6, 1070)  INS.Al 
1272  FORMAT  (1X,A10,F15.5) 

DO  4  1=1,35 

IF ( INSTR(I) .EQ.INS)  GO  TO  6 

4  CONTINUE 
STOP  4 

I  =  1  2  3  4  5  6  7  8  9  10  11  12  13  14  15  16  17  18 

6  GO  TO(ll, 11, 11, 20, 20, 20, 30, 30,33, 30, 40, 40, 50, 50. 50,50,60,60 
21  22  23  24  25  26  27  28  29  30  31  32  33  34  35 


1  70, 7 0,83,80, 81, 82, 82, 82, 83, 84, 84, 86, 86, 86, 86), I 


INPUT  THE  COLLAPSE  FIELD  HC.  KM AG  =  1 

ll' HC=A1 

NMAG=  NMAG  +  1 
KMAG(NMAG)=1 

IF(NLEN.GE.l.AND.NMAG.GE.l.ANO.NLEN*-NMAG.E2,3)  GO  TO  90 
GO  TO  1 

IN°UT  THE  STRIPOUT  FIELD  HS.  KMAG  =  2 

20  HS  =  A 1 

NN AG=NMAG*1 
KMAG ( NMAG )  =  2 

IC{NLEN.GE.1.AN0«NMAG,GE«1.ANQ.NLFN*NMAG.E3.3)  GO  TO  90 
GO  TO  1 


BUBCH1 1 5 
BUBCH1 1 6 
BU0CH117 
BUBCH1 1 8 
BUBCH1 19 
BUBCH120 
BUBCH1 2 1 
BU8CH1 22 
BU3CH123 
BUOCH1 24 
BUBCH125 
3UBCH1 2  6 
8U3CH127 
BUBCH128 
BU9CH129 
BU3CH13C 
BU9CH1 31 
BUBCH1 3  2 
B’JBCHl  3  3 
3U8CH1 34 
BUBCM135 
BUBCH1 3  6 
BUBCH1 3  7 
B'JBCHl  38 
BUBCH1 39 
BUBCH14C 
8  U  3  C  H 1 4  1 
BUBCH1 4  2 
BU3CH14  3 
3UBCH1 4  4 
B'JBCHl  4  5 
BU3CH146 
BUBCH14  7 
BU9CH1 4  8 
19  20  3J3CH149 
BUBCH15C 
60,63 , 3U3CH15 1 
BUBCH152 
BUBCH15 3 
BUBCH15  4 
3U3CH155 
BUBCH156 
BUBCH15  7 
BUBCH158 
BU9CH1 5  9 
9U3CH16G 
BUBCH16  1 
BUBCH16  2 
BUBCH16  3 
BU3CH1 6  4 
BJBCH165 
BUBCH166 
BUBCH167 
BUBCH168 
8URCH169 
3UBCH17  0 
3UBCH1 7 1 


36 


oo  ooo  ooo  o  o  o  oot)  ooo  oou  non 


AFAL-TR-79-11 28 


INPUT  THE  THICKNESS  H.  KEEN  *  1 

30  H»A1 

NLEN=NLEN»1 

KLEN|NLEN>*1 

IF<NLEN.GE.1.AM0.NMAG.GE.1.AND.NLEN»NMAG.E3.3>  GO  TO  90 
GO  TO  1 

INPUT  THE  SUSCEPTIBILITY  TIMES  THICKNESS  XuM.  KLEN  =  2 

40  X0H=A1 

NLEN=NLEN»1 

KLEN<NLEN>=2 

IF1NLEN.GE.1.AND.NMAG.GE.1.AND.NLEN+NMAG.E3.3)  GO  TO  93 
GO  TO  1 

INPUT  THE  STRIPEWI OTH  SW.  KLEN  =  3 

50  S4  -  A  1 

NLEN=NL£N»1 
KLEN(NLEN) =  3 

IFtNLEN.GE.l.AN|D.NMAG.GE.l.AND.NLENtNMAG.E3.  3)  GO  TO  90 
GO  TO  1 

INPUT  THE  SATURATION  MAGNETIZATION  MS.  KMAG  =  3 

5C  MS=A1 

NMAG  =  NMAG* 1 
KMAG  (NMAG)  =  3 

IF (NL EN.GE. 1 .ANO.NMAG.GE .1 .A ND.NLFN«-N.MAG. E3. 3 »  GO  TO  90 
GO  TO  1 

INPUT  THE  CHARACTERISTIC  LENGTH  L.  KLEM  =  4 

7  0  L  =  A 1 

NLEN=NLEN*1 
KLEN ( NLEN ) =4 

IF(NLEN.GE.1.A'JO.NMAG.G£.1.ANO.NLEN*NMAG.E3.3I  GO  TO  90 
GO  TO  1 

INPUT  THE  SUSCEPTIBILITY  XC.  KMAG  =  4 

30  XC  =  A 1 
NMAG=NMAG*1 
KMAG (NMAG ) =4 

IF  (NLEN.GE.l.AHD.NMAG.GE.l  .AN0.NLEN«-NMAG.E3.3>  GO  TO  93 
GO  TO  1 

INPUT  THE  QUALITY  FACTOR  Q  KIPM  =  1 

31  Q  =  A 1 
KIPM=1 

IF<MS.N£.0. I  GO  TO  14C 
GO  TO  1 

INPUT  THE  ANISOTROPY  FIELD  HK  KIPM  =  2 


BUBCH172 
BU9CH1 7 3 
BUBCH1 7  4 
BUBCH175 
8UBCH1 76 
BU8CH1 7  7 
BUBCH1 7  8 
BUBCH1 7  9 
8U9CH130 
BUBCH1 81 
BUBCH182 
BUBCH183 
BUBCH1 8  4 
BU3CH185 
BUBCH186 
BU3CH1 8  7 
BUBCH168 
BUBCH189 
BUBCH190 
3UBCH1 9 1 
BUBCH192 
BUBCH193 
BUBCH194 
BU9CH195 
3UBCH196 
BUBCH197 
BUOCH198 
BU3CH1 99 
BU9CH2i C 
BUBCH2  u 1 
BU3CH232 
BU9CH2C  3 
BU3CH2  C  4 
BUBCH2:5 
BU1CM20  6 
8U3CH20 7 
BU8CH?  j  8 
BUBCH2 l 9 
BU8CH2  1 3 
8UBCH21 1 
0UBCH212 
BUBC  H2 1 3 
BUBCH21 4 
BUBCH215 
8UBCH216 
BUBCH217 
BJ3CH213 
BUBCH21 9 
9UBCH22  £ 
BUBCH221 
BU9CH222 
3UBCH22  3 
BUBCM224 
BJBCH225 
BUBCH226 
BUBCH22  7 
BJBCH228 
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C 

82  H<= A1 
KI PM=2 

IF(MS.N£.C.)  GO  TO  130 
GO  TO  1 

INPUT  THE  ANISOTROPY  ENERGY  OENSITY  KU  <IPM  =  3 

83  KU  =  A  1 
KIPM=  3 

IF (MS.NE.O.  )  GO  TO  125 
GO  TO  1 

INPUT  MAGNETIZATION  DETERMINED  FROM  COLLAPSE  FIELD 
USING  THE  CALLEN  ANO  JOSEPHS  APPROXIMATION  KHAG  =  5 

84  C JM= A1 
NMAG=NMA6*1 
KMAG (NMA  G ) =5 

IF(NLEN.GE.l.ANO.NMAG.G£.l . A ND .NLEN*NM AG . EQ. 3 )  GO  TO  90 
GO  TO  1 

INPUT  AN  IDENTIFYING  LA3EL. 

85  BACKSPACE  5 
REA0(5,961>  LABEL 
GO  TO  1 

INPUT  THE  FIELD  FOR  H ALF- SATUR A T I  ON  MAGNETIZATION  KMAG  =  o 

85  HI 2=A1 

NMAG=NMAG*1 
KMAG ( NMAG ) =6 

IF(NLEN.GE.1.AND.NMAG.G£.1.AN0.NLEN«-NMAG.E3.3)  GO  TO  90 
GO  TO  1 

90  IF(KMAG(1). EQ.l. AND. KMAG(2) .EO.O)  GO  TO  101 

IFJ (KMAG (1) .EQ. 1. ANO. KMAG (2) .EQ.2) .0®. (KMAG (1) .EQ.2 .AND. KMAG (2) 
1  .EQ.l))  GO  TO  99 

IF  HKMAG (1) .EQ.  1. AND. KMAG (2) .EQ.3> .OR. (KMAG (1 ) .EQ.3. AND. KMAG (2) 
1  .EQ.l))  GO  TO  98 

IF ( (KMAG (1) .EQ. 1 . A  NO . KMAG ( 2 )  . EQ.4) .OR.  (KMAG (1 ) . EQ . 4 . A ND . KMA G ( 2) 
1  .EQ.l))  GO  TO  97 

IF(KMAGtl). EQ.2. ANO. KMAG(2) .EQ.0)  GO  TO  1C1 

IF ( (KMAG (II .EQ.2. AND. KMAG (2) .EQ.3) .OR. (KMAG ( 1 ) . EQ . 3 . AND . KMAG ( 2) 
1  . EQ • 2  I )  GO  TO  95 

IF ( (KMAG (1) . EQ. 2. AND. KMAG (2) . FQ.4 ) .OR. (KMAG (1 ) . EQ.4 .AND. KMAG (21 
1  .£Q,2>>  GO  TO  94 

IF(KKAGd). EQ.3. AND. KMAG(2I. EQ.l)  GO  TO  1G1 

IF ( (KNAG (II .EQ.  3.  AND. KMAG (2)  .EQ.4) .OR.  (KMAG (1 )  .  EQ.  4 . AND . KMAG ( 21 
1  .EQ.3))  GO  TO  92 

IF(KMAG(1). EQ.4. AND. K MAG ( 2 >  .  EQ . J  )  GO  TO  111 
IF (KMAGd ) .EQ. 5. AND.  KMAG(2 I  .  EQ. J)  GO  TO  103 
IF (KMAG(1 ) .EQ.o. ANO. KMAG(2)  .  EQ. J)  GO  TO  104 
STOP  "90 ” 

XO  ANO  MS  MERE  INPUT.  KMAG  =  3  ANO  4 


BUBCH229 

BUBCH23  u 

BUBCH231 

BU3CH232 

BUBCH233 

BUBCH234 

BUBCH235 

BUBCH236 

BUBCH23Z 

BU3CH238 

BUBCH?  39 

8UBC  H2  4  j 

BUBCH241 

BUBCH242 

BUBCH  ?43 

BUBCH2>*4 

BUBCH245 

BUDCH24G 

BUBCH?  4  7 

BU8CH24  8 

BUBCH?  4  9 

BU8CH250 

3UBCH?51 

BUBCH252 

8UBCH75  3 

8UBCH254 

3UBCH255 

BUBCH25  6 

BUBCH25  7 

BUBCH25  8 

3U3CH259 

BUBCH260 

BUBCH261 

BU8CH262 

BU3CH26  3 

BUBCH264 

BUBCH2G5 

8U8CH268 

BUDCH76 7 

BUBCH26. 8 

3UBCH2G9 

BUBCH?  7 5 

3UBCH271 

BUBCH2  7 Z 

3UBCH2Z3 

BUBCH?  7  4 

BUBCH275 

BUBCH?  7  5 

BUBCH?  7 7 

BUBCH2  7  8 

BUBCH279 

BJBCH28C 

BU3CH28 1 

BUBCH28  2 

BUBCH28  3 

BUBCH2tt4 

BUBCH285 
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92  COLL  LAGINT(X3,CHIO,OI,HC3YMS,HCOVMS,43> 

CALL  LAGINT(XO,CHIO,DI,HS9YMS,HSOVMS,40) 

CALL  LAGINT(XC,CHIO,OI,LOVH,LOH,43> 

CALL  LAGINT(XG,CHIO,OI,HOVH,HOH,4C  ) 

CALL  LAGINT(LOrt,LOVH,OI , HH ALF , H120 VM, 4 0 ) 

HC=HCOVMS*MS 

HS=HSOVMS*MS 

GO  TO  102 

X0  AND  HS  HERE  IN“UT.  KM  AG  =  2  AND  4 

94  CALL  LAGINT(XC,CHI0,0I,HSBYMS,HSOVMS,40> 

CALL  LAGINT(Xj,CHIG , 0 I , HCBYMS , HCO VMS, 4 0 ) 

CALL  LAGINT(XO,CHIO,DI,LOVH,LOH,4G) 

CALL  LAGINT(XC,CHIO,QI,WOVH,HOH,4G ) 

CALL  LAGlNT(LOH,LOVH,OI,HHALF,H12OVM,40) 
M3=HS/HSOVMS 

HC=HCOVMS*MS 
GO  TO  1 02 

HS  AND  MS  HERE  INPUT.  KM AG  =  2  ANO  3 

95  HSOVMS=HS/MS 

CALL  LAGINT (HSO VMS , HS 8YMS , 01 , HC8Y MS , HC OVMS, 4 u ) 
CALL  LAGINT(HSOVMS,HS9YMS,OI,LOVH,LOH,4C> 

CALL  LAGINT (HSO VMS, HSBYMS, 01 , WOVH, HOH, 40 > 

CALL  LAGINT (HS0VMS,HS8YMS, DI ,  C  H I  j , X . ,40) 

CALL  LAGINT (LOH.LOVH, DI.HH ALF, H120VM,4G> 

HC=HCOVMS*MS 

GO  TO  10  2 

X0  ANO  HC  HERE  INPUT.  KM  AG  =  1  ANO  4 

97  CALL  LAGINT(XO,CHI3,OI*HC8YMS,HCOVMS»40) 

CALL  LAGINT  (XO.CHIO  ,OI,HSBYMS,HSOVMS,43) 

CALL  LAGINT(XC,CHI3,OI,LOVH,LOH,4C) 

CALL  LAGINT ( XC , CHIC ,D I, WOVH, WOH, 4C ) 

CALL  LAGINT (LOH.LOVH ,DI .HHALF , HI  20 VM,40) 
MS=HC/HCOVMS 

HS=HSOVMS*MS 
GO  TO  102 

HC  AND  MS  HERE  INPUT.  KM  A  G  =  1  ANO  3 

98  HCOVMS=HC/MS 

CALL  LAGINT  (HSO  VMS,  HCOYMS,  DI  ,  WOVH  ,WIH, 4G ) 

CALL  LAGINT ( HCO VMS , HC DYMS , 01 .LOVH.LOH, 40 ) 

CALL  LAGINT (HCO VMS , HC BYMS , DI , HSBY MS , HSO VM S, 4 0 ) 
CALL  LAGINT ( HSO VMS , HC DYMS , DI , CHI 0 , X 1 , 4 3 ) 

CALL  LAGINT (LOH , LOVH , DI , HHAL F , H120 VM, 4 0 ) 

HS=HSOVMS*MS 

GO  TO  102 

HC  ANO  HS  HERE  INPUT.  KM  A  G  =  1  AND  2 

99  HCOVHS=HC/HS 

CALL  LAGINT(HCOVHS*HCBYHS,OI,LOVH,LOH,40) 


BUBCH2  8  6 
BUBCH28  7 
BUBCH288 
BUBCH289 
BU9CH290 
BUBCH29 1 
BU0CH292 
BU3CH293 
8UBCH294 
BUBCH295 
3U9CH796 
BUBCH297 
BU8CH298 
BUBCH299 
BUBCH3G  0 
BUBCH3C1 
BU3CH3G2 
BUBCH3C 3 
BUBCH3G4 
BUBCH3 l  5 
BUBCH3C6 
BUBCH3C 7 
BUBCH3C 8 
BU3CH309 
8UBCHUC 
BU8CH31 1 
BU8CH312 
BUBCH31 3 
BUBCH314 
BUBCH3 1 5 
8UBCH316 
BUBCH317 
BUBCHU8 
8U3CH319 
BUBCH32C 
OUBCH321 
BU8CH322 
8UBCH323 
BU8CH324 
BU8CH325 
BUBCH  32  6 
BUBCH3  2  7 
BU8CH328 
BUBCH  329 
BU8CH333 
BUBCH  33 1 
BU8CH332 
BUBCH  3  3  3 
BUBCH334 
BU3CH335 
BUBCH336 
BUBCH  33  7 
BUBCH338 
BUBCH3  39 
BUBCH34Q 
BUBCH341 
BUBCH342 
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CALL  LAGINT(HC3VHS,HCBYHS,DI ,W0VH, WOH.401 
CALL  LAGINT(HCOVHS,HCBYHS,OI , HCBYHS , HCOVMS, 40 ) 

CALL  LAGINT (HCO VHS , HC BYHS , 01 , HSBYMS , HSOVHS, 40  > 

CALL  L AGINT (HCO VHS, HCBYHS, 01, CHI G  ,X>.  ,40) 

CALL  LAGINT(LOH,LOVH,OI,HHALF,H12OVM,40) 
NS=(HC/HCOVMS*HS/HSOVMS) *0.5 
GO  TO  1G  2 

101  IF((KLEN(1)  .EQ.l.ANO.KLEN(2>  .  E0.2)  .OR.  (KLEN(l)  .  EQ .  2  .  A  NO  .  KLEN  (  2) 
1  .EQ.ll)  GO  TO  112 

IF( (KLEN(l) .EQ. 1.AN0. KLEN ( 2)  . EQ. 3  I . OR. (KLEN(1).EQ.3.AND.KLEN(2) 
1  .EQ.ll)  GO  TO  111 

IF C (KLEN (1) .EQ.l.ANO.KLEN(2)  .EQ.4) .OR. ( KLEN ( 1 1 . EQ . 4 . ANO . KLEN ( 21 
1  .EQ.ll)  GO  TO  110 

IFUKLEN(l)  .EQ.  2.  ANO.  KLEN  (2)  .EQ.  3)  .OR.  ( KL  EN  ( 1  I  .  EQ .  3  .  ANO  .  KLEN  <  2) 
1  . EQ. 2  1 1  GO  TO  109 

IF ( (KLEN (1) .EQ. 2. AND. KLEN (2)  .EQ.4 I  .OR.  ( KLEN ( 1 1 . EQ . 4 . A  NO . KLEN ( 21 
1  • EQ. 2 1 1  GO  TO  108 

IF( (KLEN(l) .EQ. 3. ANO. KLEN (2) . EQ. 4  I  . OR. ( KLEN ( 1 1 . EQ . 4 . A  NO . KLEN ( 2) 
1  . EQ. 3 1 1  GO  TO  107 
STOP  101 

102  IF(KLEN(1). EQ.ll  GO  TO  113 
IP(KLEN(1I.EQ.2)  GO  TO  114 
IF(KLEN(1I .EQ.3I  GO  TO  115 
IF(KLEN(1) .EQ.4I  GO  TO  116 
STOP  102 

10  3  IF ( (KLEN (II .EQ.l. ANO. KLEN (21  .  EQ . 3 }  . OR.  ( KLEN ( 1 1 . EQ . 3 . ANO . KLEN ( 21 
l.lll  GO  TO  1221 
STOP  103 

104  IF( (KLEN(l) «EQ.l.ANO.KLEN(2) .EQ.3I .OR. (KLEN (II .EQ.3.AN3.KLENI2I 
1  .EQ.ll)  GO  TO  1222 
STOP  104 

SW  ANO  L  WERE  INPUT.  KLEN  =  3  ANO  4 

107  CALL  LAGINT(-SW/L,WOVH,LOVH,LOVH,LOH,4J> 

CALL  LAGINT (L OH, LOVH, OI, HCBYHS, HCOVMS, 40  I 
CALL  LAGINT(LOH,LOVH,OI,HSBYMS,HSOVMS,4G) 

CALL  LAGINT (LOH.LOVH, OI ,CHIj ,XG,4C I 

CALL  LAGINT (LOH , LOVH , OI , HHALF , H120 VN, 4 3 1 

H=L/LOH 

GO  TO  117 

L  ANO  X0H  WERE  INPUT.  KLEN  =  2  ANO  4 

108  CALL  LAGINT(-XOH/L,CHIO,LOVH,LOVH,LOH,4G) 

CALL  LAGINT(LOH,LOVH,DI,WOVH,WOH,40 I 
CALL  LAGINT (LOH.LOYH, 01, CHIj,X0,4C) 

CALL  LAGINT (LO.H,  LOVH,  OI ,  HCBYHS  , HCOVMS, 46) 

CALL  LAGINTILOH.LOVH, OI , HS3Y MS , HSO VMS, 4C I 
CALL  LAGINT (LOH.LOVH, DI,HHALF,H120VM,4J) 

H  =  X  OH/XG 
SW=W0H*H 
GO  TO  117 

XOH  ANO  SW  WERE  INPUT.  KLEN  =  2  ANO  3 
1C  9  CALL  LAGINT (-X0H/SN, OHIO , WOVH , LO VH , LOH , 40  I 


BUBCH343 
BUBCH344 
BUBCH345 
BU9CH346 
BU3CH347 
BUBCH348 
BU0CH349 
BJBCH350 
BUBCH351 
BU3CH352 
BU3CH353 
8UBCH354 
BU0CH355 
BUBCH356 
BUBCH35  7 
3UBCH358 
3UBCH359 
BU3CH360 
BUBCH361 
BUBCH362 
BJOCH363 
BU3CH364 
BUBCH365 
BUBCH366 
3U3CH3  6  7 
EQ3U3CH36* 
BUBCH369 
9U3CH3  7  0 
BUBCH371 
BU3 CH372 
BUBCH373 
BUBCH374 
9UBCH3  7  5 
BUBCH376 
BU3CH377 
BUBCH378 
3U3CH379 
BUBCH  3  8  0 
BU3CH38 1 
0U3CH382 
BU3CH383 
BUBCH384 
BU3CH385 
BU3CH38  6 
0UBCH387 

BUBCH3  8  8 
BUBCH389 
BU3CH39G 
BUBCH (91 
BU9CH392 
BU3CH393 
BU3CH394 
BUBCH395 
BUBCH396 
BUBCH 397 
BUBCH3  9  8 
BUBCH  399 


40 
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CALL  LAGINT (LOH.LOVH, 01 , WO VH , WOH, 43  ) 

CALL  LAGINT(LOH,LOVH,OI,CHIu ,XC,40) 

CALL  L AG  INTI LOH.LOVH, 01, HC3Y MS, HCO VMS, 40 
CALL  LAGINT (LOH.LOVH, 01 , HSBY MS.HSO VMS, 40) 
CALL  LAGINT (LOH.LOVH, 01, HHALF,H120VM,4 3) 
H=X0H/X0 
SW=WOH*H 
GO  TO  11 7 

H  ANO  L  HERE  INPUT.  KLEN  =  1  AND  4 

110  LOH=L/H 

CALL  LAGINT (LOH.LOVH, 01 , HSBYMS,HSOVMS, 4C I 
CALL  LAGINT(L0H,L0VH,0I,HC3YMS,HC0VMS,4(j) 
CALL  LAGINT(LOH,LOVH,OI,CHIj,X3,4l) 

CALL  L AG INT (LOH.LOVH, 01, WO VH, WOH, 43) 

CALL  LAGINT (LOH.LOVH, 01, HHALF , H120VM, 4 0 ) 
SW=WOH*H 
GO  TO  117 

H  ANO  SW  WERE  INPUT.  KLEN  =  1  AND  3 

111  W0H=SW/H 

CALL  LAGINT(WOH»WOVH, OI.HS3YMS.HSOVNS.40) 
CALL  LAGINT (WOH, WOVH , 01 , HC3Y MS , HCO VMS, 40 ) 
CALL  LAGINT (WOH, WOVH, 01, CHI j , Xu, 40 ) 

CALL  LAGINT  (W0H,W0VH,0I,L9VH,L0H,43) 

CALL  LAGINT(LOH,LOVH,OI,HHALF,Hl2OVN,40) 
L=L0H*H 
GO  TO  117 

H  ANO  XOH  WERE  INPUT.  KLEN  =  1  ANO  2 

112  XO  =  XC  H/H 

CALL  LAGINT ( XC  ,  CHIO  « D I , HS3 YMS , HSO VMS , 4 3 ) 
CALL  LAGINT  t  XC  ,  C  HIO  ,  0  I ,  HCO  YMS ,  HCO  VMS,  *»  0  ) 
CALL  LAGINT(XC,CHIC,OI,WOVH, W0H.4C ) 

CALL  LAGINT(XC,CHIO,OI,LOVH,LOH,4ij  I 
CALL  LAGINT  (LOH.LOVH,  01,  HH  ALF,  H12  0VM,  i,0> 
L=LOH*H 
SW=WOH*H 
GO  TO  117 

H  WAS  INPUT  KLEN  =  1 

113  L-LOH*H 
SW=WOH*H 
GO  TO  122 

XC H  WAS  INPUT  KLEN  =  2 

114  H=XCH/X3 
L=LOH*H 
SW=WOH*H 
GO  TO  122 

C 

C  SW  WAS  INPUT  KLEN  =  3 


BU3CH40  0 
BU3CH4C1 
BU3CH4C2 
BUBCH40  3 
BU0CH4  3  4 
BU3CH4C5 
BUBCH40  6 
BUBCH40  7 
BUBCH40 8 
BU8CH4C9 
3UBCH41Q 
BU3CH4 1 1 
BUBCH412 
BU9CH41 3 
BUBCH414 
BUBCH415 
BUBCH4 16 
BUBCH4 1 7 
BU3CH418 
8U3CH4 1 9 
BU3CH4  2  0 
3U3CH42 1 
BU3CH422 
BU3CH42  3 
BU3CH42  4 
BUBCH4  2  5 
BUBCH426 
BU3CH42  7 
3U3CH428 
8UBCH429 
BUBCH433 
BUBCH4  3 1 
BUBCH43  2 
BU3CH4  3  3 
BUBCH434 
3U3CH435 
BUBCH<«36 
BUBCH43  7 
8UBCH438 
8U8CH439 
BU9CH44  C 
BUBCH44 1 
BU3CH442 
BUBCH44? 
B'JBCH',  4  4 
8U9CH445 
BU3CH446 
9U9CH44  7 
BU3CH448 
BUBCH449 
8U9CH46C 
BU3CH451 
BUBCH462 
8U3CH46  3 
BUBCH4S4 
BUBCH46S 
BUBGH466 
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C 

115  H  =  S  H  /  H  OH 
L=LOH*H 
GO  TO  122 

L  HAS  INPUT  KLEN  =  4 

116  H=L/LOH 
SH=HOH*H 
GO  TO  122 

SELECT  THE  MAGNETIC  PARAMETERS  YET  TO  BE  CALCULATED 

117  IF(NMAG.GT.l)  STOP  117 
IF (KMAG* 1 ) . ED . 1 )  GO  TO  118 
IF (KMAG(1 > . E0.2)  GO  TO  119 
IF(KMAG(1) .EQ.3)  GO  TO  121 
IF(KMAG(1).EQ.4J  GO  TO  121 
STOP  1170 

HC  HAS  INPUT  KMAG  =  1 

118  M3=HC/HCOVMS 
HS=HSOVMS*MS 
GO  TO  122 

HS  HAS  INPUT  KMAG  =  2 

119  MS=HS/HSOVMS 
HC=HCOVMS*MS 
GO  TO  122 

MS  HAS  INPUT  KMAG  =  3 

120  HC=HCOVMS*MS 
HS=HSOVMS*MS 
GO  TO  122 

X0  HAS  INPUT  KMAG  =  4 

121  HRITE (6, 1010)  KC , H , S W , L , HC 0 VMS, HSO VMS 

CIO  FORMAT  (/"  X0  HAS  THE  ONLY  MAGNETIC  PARAMETER  INPUT"/"  X j  =  "E12 

1  ",  H  =  "E12.6",  SW  =  "£12.f>",  AND  L  =  "E12.6"."/"  ALL  THAT  CAN 

2  SAID  ABOUT  THE  OTHER  MAGNETIC  PARAMETERS  IS*"/"  HC/4.PI.MS  =  " 

3  E12.6"  AND  HS74.PI.MS  =  "E12.6"."//) 

GO  TO  10 

CALLEN  AND  JOSEPHS  MAGNETIZATION,  STRIPEH I OTH,  AND  THICKNESS 
HERE  INPUT.  CALCULATE  HC,  SET  KMAG  TO  1,  ANO  GO  TO  111. 

1221  H0H=SH/H 

CALL  LAGINT (HOH,HOVH, DI,LOVH,LOH,40 ) 
HC=(1.+0.75*LOH-SORT(3.*LOH) ) *CJM 
KMAG  1 1  )  =1 
GO  TO  111 
C 

C  LASER  SPATIAL  FILTERING  CHARACTERIZATION. 


BU3CH457 
BU8CH458 
BUBCH459 
BU3CH46 j 
BUBCH461 
BUBCH462 
BU8CH463 
BU8CH46U 
BUBCH465 
BU9CH466 
BU3CH467 
BU8CH468 
BUBCH469 
BUBCH470 
BUBCH4  7 1 
BUBCH472 
BUBCH473 
BU3CH474 
3UBCH475 
BU3CH4  7 6 
BUBCH4  77 
BUBCH478 
BUBCH4  7  9 
8U9CH480 
BU3CH4  8 1 
BUBCH482 
3U9CH48  3 
BUBCH484 
BU3CH48S 
BUBCH486 
BUBCH487 
BUBCH488 
BUBCH4  8  9 
BU3CH491 
BU3CH49 1 
3UBCH492 
BUBCH493 
BUBCH494 
3U3CH495 
8U3CH49  6 
8U0GH497 
?>  BU3CH498 
3E0U3CH499 
BU3CH3  o  0 
8UBCH5  u 1 
BUBCHB  C 2 
BU3CH5C  3 
BU3CHBC4 
BU3CH50  5 
BU3CH506 
BUBCHSu  7 
BUBCH5C8 
BU3CH5C  9 
8U9CH5 1 0 
BU9CH5 11 
BUBCH5 1 2 
BU3CH5 1 3 
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FIELD  (H12)  FOP  HALF-SATURATION  MAGNETIZATION, 

THICKNESS,  ANO  ZERO-MAGNETIZATION  STRIPEWIOTH  HERE  INPUT. 

CALLEN  AND  JOSEPHS  APPROXIMATE  ANALYSIS  IS  ASSUMED  AND 
CORRECTED  BY  THIELE  ANALYSIS  USING  OERIVEO  BUBBLE  COLLAPSE 
FIELD 


1222  W3H=SW/H 

CALL  LAGINT (WOM, WOVH, 01,  HHALF,H120VM,4j) 

CJM=H12/H120VM 
GO  TO  1221 
122  CONTINUE 

H12  =  H120I/M*MS 
WALL=MS*MS*L/PI4 
WALL=WALL*1 .E-4 

WRITE (6, 1C  00)  LABEL, H,L, X0,SW,MS,HS,HC  .WALL, HI  2 
lOuC  FORMAT  (/IX, 8Alj/”  H  =  "E12.6",  L  =  "E12.6",  CHIJ  =  "E12.6",  SW  = 

1  El 2 . 6/"  4. PI. MS  =  "E12.6",  HS  =  "E12.6",  HC  =  "E12.6",  BULK  WALL  BU3CH531 
2ENERGY  =  "E12.5/”  Hl/2  ( FI  EL  0  FOR  H ALF-SA T J RA T ION  MAGNETIZATION)  =3UBCH532 


BU3CH51 4 
BU8CH515 
BU3CH516 
BUBCH51 7 
3U8CH51 8 
BUBCH519 
BUBCH52  0 
BU8CH521 
BU3CH52? 
BUBCH52  3 
BUBCH524 
BUBCH525 
8U3CH525 
BU8CH52  7 
BUBCH528 
3UBCH529 
3U3CH530 


3  “E 1 2 . 6//  ) 

9JBCH533 

CALL  STA9LT (LOW, H, MS) 

BU3CH534 

IF(KIPM.NE.O)  GO  TO  124 

BUBCH535 

READ ( 5 ,960 )  INS 

BUBCH536 

IF(E0F(5)I  100,123 

BU3CH53  7 

123 

BACKSPACE  5 

BU3CH5  3  8 

IF(INS.EQ."Q".OR.INS.EQ."HU”.OR. 

INS.£Q."HA“ 

.OR. INS. EQ." 

HK”. OR .  BUBCH5  39 

1  I NS . EQ. “KU“ )  GO  TO  10 

BU3CH54 j 

GO  TO  9 

BU3CH541 

124 

CONTINUE 

BUBCH542 

WRITE (6, 1090 ) 

BU9CH54  3 

1090 

FORMAT (/10X“IN-PLANE  PARAMETERS 

(ACCORDING 

TO  DEBONTE) 

FOLLOW!"// ) BJ3CH544 

IF(KIPM.EQ.l)  GO  TO  140 

BU3CH545 

IF(KIPM.EQ.2)  GO  TO  130 

BU3CH546 
BUBCH54  7 

ANISOTROPY  ENERGY  DENSITY  KU 

WAS  INPUT 

KIPM  =  3 

BU3CH548 

BUBCH549 

125 

CONTINUE 

BU3CH55C 

HK=2. *KU*PI4/MS 

BUBCH55 1 

Q=HK/MS 

BUBCH552 

GO  TO  150 

3U3CH55  3 
BU3CH554 

ANISOTROPY  FIELD  HK  WAS  INPUT 

KIPM  = 

2 

BUBCH555 

BU3CH55b 

130 

Q=HK/MS 

3UBCH557 

KU=j.5*MS*HK/PI4 

BU3CH55  8 

GO  TO  150 

BUBCH559 

BU3CH56C 

DUALITY  FACTOR  D  WAS  INPUT 

KIPM  =  1 

3UBCH561 

3UBCH562 

140 

MK=Q*MS 

BU3CH56  3 

KU=0.5*MS*HK/PI4 

BUBCH564 

BUBCH565 

COMPUTE  THE  REMAINING  INPLANE 

PARAMETERS 

0UBCH566 

BUBCH567 

150 

EXC0NS=WALL*WALL/(16.*K'J) 

BUBCH568 

COASHl=Q/LOH 

BUBCH589 

CALL  LAGRGE(COASHl,PHll,Yl) 

BUBCH5  7  0 
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CAUL  LAGRGEJG0ASH1»EQASH,Y2) 

CALL  LAGRGE JC04SH1, RATIO, Y3» 

ALENTH=JPI4*0.25/Y3)*SQRTJEXCONS/KU) *10000. 

YU  =  Y1*720./PI4 
Y2=Y2*2.C*SQRTIKU*EXCONS) 

WRITE J6.1C83 )  fll.Yl ,Y2, AL £N TH, Q, KU , HK , E XCONS 
1080  FORMAT!////"  MAXIMUM  TWIST  ANGLE”13! ". ")"="F8. 4* 3X "DEGREES “3X« 

1  “!~F5.2"  RAO.)”/"  DOHA  I N  HALL  ENERGY"14 ! " . ") ”="E 14 . 6 . 3 X, 

2  3X“ERGS/SQ.CM.“/“  HALL  THICKNESS "18!",”)“= "£14. 6, 3X 

3  “MICRONS"/”  3  =  "Flu. 3",  KU  =  "E12.6",  H<  =  “E12.6",  AND  EXCHA 
4E  CONSTANT  =  ”212.6) 

MS  =  0  • 

GO  TO  IQ 
100  STOP 
END 

SUBROUTINE  LAGI NT ( X , A RG1 , A RG2 , tf AL . Y , NO IM) 

THIS  13  THE  SUBROUTINE  WHICH  PERFORMS  THE  4-POINTS 
LAGRANGES  I  NT ERPRE T A T ION  FOR  THE  OATA 

01  MENS  ION  ARGK56)  ,  ARG2(56)  ,VAL!56),RR!56) 

DO  10  J=  1 »  5  6 
10  RR!J)=0. 

IFJX.LT.j.)  GO  TO  6 
OO  1  I=1,N0IM 

IF (ARGl(I) . EO . 0 . . A  NO  .  X  .  NE . 0.  )  GO  TO  1 
IF(ARG2(I) .Lt.l.E-6)  ARG2 < I ) = 1 .  E- 6 
IF!ARG1!I)/ARG2!I)-X)  1,2,3 

1  CONTINUE 

2  Y  =  VAL ( I) 

RETURN 

3  IPIV=I 


BU8CH5  7 1 
BUBCH572 
3UBCH573 
3UBCH574 
BUBCH575 
BU3CH576 
3UBCH577 
BUBCH578 
3U3CH579 
NG3UBCH580 
BU3CHS  8 1 
BUBCH58  2 
BU3CH583 
BU3CH584 
BUBCH585 
LAGNT  1 
LAGNT  2 
LAGNT  3 
LAGNT  4 
LAGNT  5 
LAGNT  6 
LAGNT  7 
LAGNT  8 
LAGNT  9 
LAGNT  1C 
LAGNT  11 
LAGNT  12 
LAGNT  13 
LAGNT  14 
LAGNT  15 
LAGNT  16 
LAGNT  17 


IF (IPIV.EQ.2)  M I  N  =  1 
IFIIPIV.GT.2)  M I N= IPI V-2 
MAX=IPIV*1 
IFJMIN.LT. 1)  MI N  =  1 
IF(MIN.GT.NOIM)  MIN=NDIM 
IF(MAX.GT.NOIM)  MAX=NOIM 
IFJMAX.LT. MIN)  MAX  =  M  I  N 
DO  4  N=MI N, MAX 
RRJN) =1. 

00  4  I=MIN,MAX 
IFJN.EQ.II  GO  TO  4 

IF  JABS  JARG1  <N)  /  ARG2  JN)  -  ARG1  J  U/ARG2  JI 
RR  JN) =  RR JN) *  JX-ARGl 1 1 )/ARG2 ( I ) ) / J  A  RG1 

4  CONTINUE 
Y  =  0  • 

00  5  I=MIN,MAX 

5  Y=Y*RRJI)*VALJI) 

RETURN 

6  X  =  A  f)S  JX) 

00  7  1=1, NO I M 

IF  (ARGKI  )  .ECl.j.  .ANO.X.NE.  3.  )  GO  TO  7 
IF JARG2JI ) .LT.l .E-6)  ARG2  J I ) = 1 . E- 6 
IP  JARG1JD/ARG2  (II-XI  3,2,7 

7  CONTINUE 

STOP"  SUBROUTINE  LAGINT  FAILED" 


LAGNT  18 
LAGNT  19 
LAGNT  2 G 
LAGNT  21 
LAGNT  22 
LAGNT  23 
LAGNT  24 
LAGNT  25 
LAGNT  26 
LAGNT  27 
LAGNT  28 

).LE.l.E-6)  GO  TO  4  LAGNT  29 

N)/ARG2(N)-ARG1  (D/ARG2JI))  LAGNT  3C 

LAGNT  31 
LAGNT  32 
LAGNT  33 
LAGNT  34 
LAGNT  35 
LAGNT  36 
LAGNT  37 
LAGNT  38 
LAGNT  39 
LAGNT  40 
LAGNT  41 
LAGNT  42 


44 
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SUBROUTINE  ST A3L T U OH , 0, SA T> 

THIS  IS  THE  SUBROUTINE  WHICH  GIVES  THE  OPERATING  POINTS 
FOR  THE  BUBBLE  DOMAINS  AS  WELL  AS  THE  COLLAPSE  ANO  THE 
RUNOUT  OIAMETERS  ANO  FIELDS 

DINENSION  01(56) 

DIMENSION  A  (56) ,F<56>  ,S0 ( 56 1 , S2 ( 5 6 >  ,  DIM ( 15) , H APP ( 15 ) 

REAL  LOH 
DATA  01/06*1./ 


LAGNT  43 
STA9L  1 
STABL  2 
STABL  3 
STABL  4 
STABL  5 
STABL  6 
STABL  7 
STABL  8 
STABL  9 
STABL  1C 


DATA  A /G. CO, 0.10,0. 20, 0.30, 0.40,0. 50, ii. 60,0. 7^,3. 80, L. 90, 1.0C.1.1CSTABL  11 

1  tl.2tl.3.1.4*1.5tl.6tl.7yi.8tl* 9.2. *2. 1.2. 2. 2.3.2. 4. 2. 5.2. 6.  STABL  12 

2  2.7t2.8t2.9t3.t3.2t3.4*3.6t3.8t4.t4.2t4.4t4.6t4.8t5.t5.2t5.4t  iTABL  13 

3  5.6»5.8, 6. ,6.2,6. 4,6.6, 6.8,7. ,7. 2, 7, 4, 7. 6, 7. 8, 8,/  iT A 8L  14 

3  DAT  A  F/u.t. 0939, .1765, .2493, .  3137, .  3708, .4216, .4672, .5083.  STABL  15 

1  .5455, .5794, .6104, .6393, .6655, .6901, .7130, .7345, .7547, .7737,  STABL  16 

2  .7917,. 8C87, .8249, . 8404, . 3551,. 8692, . 8827, .8956, .9081, .9200 ,  STABL  17 

3  .9316, .9427, .9639, .9837, 1.3024,1. 0 ?C 1 , 1 . u3 68 , 1.0526,  1.  06  73,  STABL  18 

4  1.0822,1.0960,1.1092,1.1219,1.1341,1.1453,1.1572,1.1681,1.1787,  STABL  19 

5  1.1889,1.1988,1.2084,1.2177,1.2268,1.2356,1.2442,1.2525,1.26.6/  STABL  20 
DATA  S0/C ... 0u59, . 0215, . 0442 , .0 716, . 1017, .1 332, . 1648, .1960, .2262  STABL  21 

C  ,.2552  STABL  22 

1  ,.2829, . 3093, .3343, .3579, .38C4, . 4016, .4219 , .4410 ,. 4592, .4765,  STABL  23 

2  . 49 31, .50 89, .52 40, .5385, .55 24, .5657, .5786, .5939, .6 028, .6143, .6 36 2STABL  24 

3  ,  .6566,  .6759, .6940, .  7112,  .  7275,  .  74’.-,  .  7573, .  7720,.  7855, .  7935,  STABL  25 

4  . 8109,. 8229, . 9345, . 3456, . 8564, . 8668 ,. 8769, . 9866, . 9961, . 9053,  STABL  26 

5  .9142, .9229, .9314, .9396/  STABL  27 

DAM  S2/0. ,.0007, .0028, .0063, .0111, .3172, .0243, .0323,. 3411, .C5C5,  STABL  28 

1  0.0603, .0705, .0809, .0914, .1020, .1126, .1231, .1336, .1439,. 1541,  STABL  29 

2  .1642, . 1741, . 1338, .1933, .2027,. 2119, .22j9,. 2297, .2838, .2458,  STABL  30 

3  .2551 ,. 2712, .2866, . 3051 ,. 3158, . 3295, . 3428, . 3556, . 3679, .3799,  STABL  31 

4  . 39 14,. 40 26, .4134, .4239, .*341, . 4439, . 4535, .4629, .472C , .4808, .48 94STA9L  32 


5  ,  .4978, .5060, .  5140, .5218, .  5295/ 

0=0*1. E-4 
HL=l./LOH 
HHL=1 . /HL 

CALL  LAGINT (HH_ ,SC ,OI ,A , X0 ,56) 

CALL  LAGINT (HHw,S2,DI,A,X2,56) 

STEP= ( X2-X0) /ID . 

WRITE (6, 1C) 

FORMAT  (IX, 52 ("•“),/»“  BUBBLE  DIAMETER  (MICRONS) ”,1CX, 
1  “APPLIED  FIELO  (OE) ",/, IX, 52 ("--) ) 

OIN ( 1 1 «Xfl  *D 
ASPECT =X0 
OO  20  1=1,11 

CALL  LAGINT (AS5tCT,A,ni,F,XF,56) 

HAPP( I MSAT ‘(YF-HHL) /ASPECT 
Oim  =  DIM(I)  *lt000. 

HAPPI =MAPP( I ) 

IF(I.EQ.l.OR.I.EQ.ll)  GO  TO  35 
WRITE  (6,30)  OHI, HAPPI 
FORMAT  (11X,F6.2,24X,F6.2) 

ASPECT =AS PE CT» STEP 

0IM(I*1)=ASPECT*0 

OIM(l)=OIH(l)*tCiC9. 

OIM(11)=DIM(11)*1C003. 


STABL  33 
STABL  34 
STABL  35 
STABL  36 
STABL  37 
STABL  38 
STABL  39 
STABL  40 
STABL  41 
STABL  42 
STABL  43 
STABL  44 
STABL  45 
STABL  46 
STABL  47 
STABL  48 
STABL  49 
STAHL  50 
STABL  51 
STABL  5? 
STABL  53 
STABL  54 
STABL  55 
STABL  56 
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WRITE (6,40)  DIN ( 1 ) »  HAPP ( 1 ) , DIM  (1 1 ) .HAPP(ll) 

40  FORMAT  (1X,52("-")  *//,•"  BUTTLE  COLLAPSE  DIAMETER  =  ",Ff>.2,2X, 

1  “MICR ON S",/«2X," BUBBLE  COLLAPSE  FIELD  =  "F6.2,"  OERSTEOS" 

2  “  BUBBLE  RUNOUT  DIAMETER  =  ",F6.2,"  MICRONS",/, 

3“  BUBBLE  RUNOUT  FIELD  =  ",F6.2,“  OERSTEOS"5  (/ ) ) 

D=0*1. E4 
80  RETURN 
ENO 

SUBROUTINE  L AGRGE ( X , VAL , Y ) 

THIS  IS  THE  SUBROUTINE  WHICH  PERFORMS  THE  4-POlNTS 
LAGRANSES  INTERPOLATION  FOR  THE  DATA  OF  OFBONTE.. 

DIMENSION  ARG ( 30 ) ,VAL (120) ,RR(3u ) 

COMMON  /AA/QU 

IN  SUBROUTINE  LAGRGE ,  "ARG"  IS  THE  DEBONTE  FUNCTION  C' 

DATA  ARG/. 1, ,13, .2, .3, ,4, .5, .6, .7, .8, ,9,1., 1,6, 2,, 3,, 4. ,S. ,6. 
1  ,8.,9.,lG,,13,,20,,30,,40.,3t.,60.,7u.,8G,,9u,/ 

DO  10  J=  1 , 30 
10  RR(J)=0. 

DO  1  K=1 , 30 
I  =K 

IF(ARG(I)-X) 1,2,3 

1  CONTINUE 

2  IWW=1 
GO  TO  6 

3  IWW=2 

6  IF(QU.LT. 1.2ft)  J=0 

IF( QU. G£. 1.28. AND. QU.LE. 1.65)  J=30 
IF (QU.GT .1.65.AN0.0U.LE. 3. 0)  J=60 

IF (OU. GT. 3. )  J=90 
IC(IWW.EQ.2)  GO  TO  7 
V  =  VAL  ( I* J ) 

GO  TO  8 

7  IFJI.LT.3)  GO  TO  74 
MIN=I -2 
M AX  =  I  +  1 

IF (MAX  .GT. 30 )  30  TO  76 
GO  TO  73 
74  MI N= 1 
MAX  =  4 
GO  TO  78 
7b  M A  X  =  3  0 
MIN=27 

78  CONTINUE 

DO  4  N=MIN,MAX 
RR(N) =1, 

00  4  I=MIN,MAX 
IF(N.EQ.I)  GO  TO  4 

I F ( ABS(ARG(N) -ARG( 1)1 .LE.O . 00001)  30  TO  4 
RR(N)=RR(N)*((X-ARG(I))/(ARG(N)-ARGtI))) 

4  CONTINUE 
T  =  0.G 

DO  5  I =MI N , M A X 
M  =  J»I 


STAOL  57 
STABL  53 
/ .STABL  59 
STABL  60 
STABL  61 
STABL  62 
STABL  63 
STABL  64 
LAGRG  1 
LAGRG  2 
LAGRG'  3 
LAGRG  4 
LAGRG  5 
LAGRG  6 
LAGRG  7 
LAGRG  8 
LAGRG  9 
LAGRG  10 
7. LAGRG  11 
LAGRG  12 
LASRG  13 
LAGRG  14 
LAGRG  15 
LAGRG  16 
LAGRG  17 
LAGRG  1ft 
LAGRG  19 
LAGRG  20 
LAGRG  21 
LAGRG  22 
LAGRG  23 
LAGRG  24 
LAGRG  25 
LAGRG  26 
LAGRG  2’ 
LAGRG  28 
LAGRG  29 
LAGRG  30 
LAGRG  31 
LAGRG  32 
LAGRG  33 
LAGRG  34 
LAGRG  35 
LAGRG  36 
LAGRG  37 
LAGRG  38 
LAGRG  39 
LAGRG  40 
LAGRG  41 
LAGRG  42 
LAGRG  43 
LAGRG  44 
LAGRG  45 
LAGRG  46 
LAGRG  47 
LAGRG  <<8 
LAGRG  49 
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LAGRG  50 
LAGRG  51 
LAGRG  5? 
LAGRG  53 
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